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SUMMARY 


TFaNS is the Tone Fan Noise Design/Prediction System developed by Pratt & Whitney under contract 
to NASA Lewis. The purpose of this system is to predict tone noise emanating from a fan stage 
including the effects of reflection and transmission by the rotor and stator and by the duct inlet and 
nozzle. These effects have been added to an existing annular duct/isolated stator noise prediction 
capability. 

The underlying concept for the system was presented in Reference 1 with application to cascades in 
2- dimensional channels. TFaNS extends this to annular geometry with duct terminations and 
radiation to the far -field via the following scheme: “Acoustic elements” (e.g. the inlet, the rotor, the 
stator, and the nozzle) are first analyzed in isolation to determine their modal reflection and 
transmission coefficients (including frequency scattering in the case of the rotor). Then the elements 
are coupled as a linear system via the duct eigenmodes at “interface planes” separating the elements. 
The linear system is solved to find a “state vector” of mode amplitudes at the interface planes. The 
“state vector” then is used to compute upstream and downstream modal sound powers and the sound 
pressure directivities in the outside field. 

TFaNS consists of: 

• The codes that compute the acoustic properties (reflection and transmission coefficients) of the 
various elements and writes them to Acoustic Properties Files, 

• CUP3D: Fan Noise Coupling Code that reads these files, solves the coupling equations, and 
outputs the desired noise predictions, 

• AWAKEN: CFD/Measured Wake Postprocessor which reformats CFD wake predictions and/or 
measured wake data so they can be used by the system to predict noise. 

Acoustic properties can be computed from a variety of codes other than those presently in TFaNS. 
For example, rotor and stator reflection and transmission are currently being computed from a 
classical, uniform axial flow model with solid body swirl. However, more sophisticated codes (e.g. 
linearized Euler codes) are anticipated for future application. CUP3D has been written in a general 
form so that acoustic properties from a variety of codes can be handled. To accomplish this, a standard 
file format for acoustic properties has been developed and must be followed. 

This document provides technical background for TFaNS including the organization of the system and 
CUP3D tec hni cal documentation. This document also provides information for code developers who 
must write Acoustic Properties Files in the CUP3D format. 

As described on page 3, technical documentation of the other TFaNS codes may be found in 
References 2 to 7. The TFaNS user’s manual may be found in Reference 8. Evaluation of the system 
may be found in Reference 9. 
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1. INTRODUCTION 


The TFaNS coupling and noise prediction scheme is explained conceptually with reference to 
Figure 1. 


Far-Field 

Microphone 



1 2 3 Interface Planes 

A 1 A 2 A 3 State Vector Elements 

B 1 B 2 B 3 Source Vector Elements 


Figure 1: TFaNS Coupling and Noise Prediction Scheme 

The fan stage is divided by three interface planes into four acoustic elements (inlet, rotor, stator, and 
nozzle). State vector, ,4, is in three sections. A 1 , A 2 , and A 3 , whose elements are modal amplitudes that 
will be defined in a later section. Source vector, B, has the same structure. B is prescribed 
(corresponding, for example, to the output of a stator due to rotor wake input in an uncoupled 
environment) and A is to be found as a function of B by solving the linear system equations. 

Coupling of the elements at the interface planes is specified in terms of the scattering matrix, S, which 
is built up from modal reflection and transmission coefficients. In condensed form, the system 
equations are represented by 

A = SA+B (1) 

which is to say that the state vector elements are the sum of the parts from scattering, SA, and directly 
from the source, B. These equations are solved formally by 

A^il-Sy'B ( 2 ) 
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Two major forms of noise output are computed from the state vector: 

• Upstream and downstream propagating sound power level are computed from 
A on a mode-by-mode basis. Power is calculated just upstream and just 
downstream of the noise source (defined in Figure 1 by interface plane 1 
upstream, and interface plane 3 downstream). 

• Outside far— field directivity is computed from the elements of A. In this case, 
far -field directivity shapes are computed by the radiation codes with unit 
amplitude input and stored as part of the acoustic properties files. These 
directivities are then multiplied by A to get the far— field sound pressure level 
directivity. 

TFaNS (Version 1.4) consists of the following computer codes: 

• CUP3D: Fan Noise Coupling Code Version 2.1 (Section 3., Appendices I to IV) 

• SOURCE3D: Rotor Wake/Stator Interaction Code Version 2.5 (Reference 2) 

• Eversman Inlet Radiation code Version 3.0 (References 3, 4, 5 and Appendix V) 

• Eversman Aft Radiation code Version 3.1 (References 4, 5, 6 and Appendix V) 

• AWAKEN: CFD/Measured Wake Postprocessor Version 1.0 (Reference 7) 

The following sections present the organization of TFaNS, the technical background behind the 
CUP3D code, and the information necessary to write an Acoustic Properties File. The first part of 
the document discusses the TFaNS organization. Then the technical documentation for CUP3D is 
presented. This includes background on the coupling system, state vector conventions, definitions of 
the reflection and transmission coefficients (scattering coefficients) in terms of the state vector, and 
the calculation of far-field directivities and mode power levels. Appendix I presents the information 
necessary to create Acoustic Properties Files. Changes made to the Eversman inlet and aft radiation 
codes to make them part of TFaNS are also presented in Appendix V. 
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2. GENERAL ORGANIZATION OF TFaNS 


The organization of TFaNS is illustrated in Figure 2. This figure identifies the codes which comprise 
the system and how they interact with each other. 

The central portion of the system is the CUP3D Fan Noise Coupling Code. This code reads Acoustic 
Properties Files which contain scattering (transmission and reflection) coefficients from other codes 
along with far— field directivity shapes and source vector information (e.g. noise from rotor 
wake/stator interaction). This information is used to form a system of linear equations which permit 
acoustic elements to reflect and transmit to each other. A System File is also input which determines 
the organization of the acoustic elements. Output from this code includes far-field directivities along 
with inlet and aft power levels. 


AWAKEN 

CFD/Measured Wake Postprocessor 

SOURCE3D Input File 
(with CFD or Measured Wakes) 

I 


Acoustic 

Wake/Turbulence 

File 


Eversman 

In let Radiation Cod e 
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Acoustic 

Properties File 
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Rotor Wake/Stator Interaction Code 

I Ro tor | 1 Stator] 

Acoustic Acoustic 

Properties File Properties File 


Eversman 
A ft Radiation Co de 
| Aft Radiation] 
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Properties File 



Figure 2: Organization of TFaNS Version 1.4 

The SOURCE3D Rotor Wake/Stator Interaction Code is a significantly extended and improved 
version of the V072 Rotor Wake/Stator Interaction Code (Reference 10, 11 and 12). It has two 
functions within TFaNS: firstly, it calculates tone noise from a rotor wake/FEGV interaction and, 
secondly, it determines the scattering coefficients for the rotor and stator then outputs them to rotor 
and stator acoustic properties files (Figure 2) for use by CUP3D. This code can either use its own 
internal semi -empirical wake model, or it can use CFD or measured wakes processed through the 
AWAKEN CFD/Measured Wake Postprocessor. 

The AWAKEN CFD/Measured Wake Postprocessor creates a SOURCE3D input file which contains 
upwash wake harmonic amplitudes calculated from CFD predictions or measured velocity data. CFD 
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or measured velocity information is obtained from the Acoustic Wake/Turbulence File which is 
generated either by a CFD code (or postprocessor) or during a engine/ng test program. 


The Eversman inlet and aft radiation codes are run if far-field directivities are required TTiese codes 
comprise three “modules”. The first module creates a finite element mesh for the calculation. The 
second module calculates the potential steady flow field. Finally a radiation module, modified to 
interface with TFaNS (see Appendix V), is used to complete the potential steady flow calculation for 
a given duct flow condition and calculate the far field radiation and scattering coefficients for a 
specified number of blade passing frequency (BPF) harmonics on a mode-by-mode basis. 
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3. CUP3D FAN NOISE COUPLING CODE 


3.1 GENERAL ORGANIZATION OF CUP3D 

The CUP3D Fan Noise Coupling Code is the final link in the TFaNS Tone Fan Noise design/prediction 
System. Output files from isolated elements are used by CUP3D to “couple” the elements thus 
accounting for reflection and transmission of acoustic and vortidty waves throughout the system. 

Inlet Rotor Stator Aft 

Acoustic Acoustic Acoustic Acoustic 



Figure 3: CUP3D File Organization 


Acoustic Elements (or Regions) 



Figure 4: CUP 3D Geometry and Terminology 
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Figure 3 is a block diagram showing the files required by CUP3D. First, codes such as an inlet 
radiation code, rotor code, stator code, and aft radiation code are run. These codes create Acoustic 
Properties Files which contain all the acoustic information about an acoustic element. The 
information may include geometry and performance, harmonic and mode numbers, scattering 
coefficients, source vector coefficients and far-field directivity shapes. Acoustic Properties Files can 
be output for any rotor, stator or radiation code. These Acoustic Properties Files and a user generated 
Systena File are then input into CUP3D. The System File contains information on how the various 
acoustic elements will be organized to form the total acoustic system. An example of an acoustic 
system is shown in Figure 4. 

Information about the terminology and geometry used by CUP3D is given in Figure 4. The present 
system can be organized with up to four acoustic elements (or regions). These elements are classified 
by any of four types: inlet radiation, rotor, stator, and aft radiation. Hie number and types of elements 
used may vary from case to case. However, if radiation elements are not used, they must be replaced 
by non-reflecting regions which act to propagate the noise away from the system. Output from 
CUP3D includes inlet and aft mode and total power levels and far— field directivities. 

Each acoustic element is bounded by interface planes. Rotor and stator acoustic eiements are 
bounded by two interface planes. Inlet and aft radiation acoustic elements are bounded by one 
interface plane and the far-field. These boundaries are non-reflecting and permit waves to be sent 
into an acoustic element one at a time. Reflection and transmission scattering coefficients can then 

be calculated. 


3.2 


STATE VECTOR AND COUPLING COEFFICIENT CONVENTIONS 


Reference 1 discusses in detail the physics of mode scattering in two dimensional flow. The purpose 
of this section is to present the equations needed to create scattering coefficients and source vector 
coefficients for the CUP3D code which handles information in three dimensions. This section is 
presented for the benefit of the source or radiation code developer who may need to create an 
Acoustic Properties File. 

“Coupling coefficient” is the collective term which refers to reflection and transmission coefficients. 
These are the matrix elements in the scattering matrix, S, that specify how any mode input to an 
acoustic element is transmitted and reflected and scatters into other modes. Since these coefficients 
are defined in terms of state vector elements, or modal amplitudes, the modal representation is 
specified first in this section. Then the coupling coefficients are defined in terms of them. 

Figure 5 shows the coordinate systems for CUP3D. Inside the engine the coordinate system is a 
cylindrical (r, <p, x) right handed system where x is positive m the downstream direction. The far held 
uses a polar coordinate system (R, <p, 8) where 6 is zero in front of the inlet. 

As in Reference 1, the flow perturbations are represented by a series of pressure and vorticity waves_ 
In general, the wave types, W, which are scattered by an acoustic element are affected by the type of 
flow through the duct. CUP3D uses a general set of modes based on whatever flow conditions are 
assumed at the interface planes (e.g. Reference 14). Though much of this report is written assuming 
a specific set of modes (given below), CUP3D is not limited to this set of modes. One exception: Sound 
power in CUP3D is presently computed using a specific set of acoustic modes. This calculation may 
be generalized in later code versions. 
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Figure 5: CUP3D Coordinate Systems 

Thus, for the purposes of the present TFaNS system, consider a constant area annular duct with a 
uniform axial flow and swirl which takes the form of solid body swirl (i.e. V 9 = rr). For this more 
limited case the wave types are: 

W = 1 Upstream going pressure wave 
W = 2 Downstream going pressure wave 

W= 3 Vorticity wave 1 (downstream only) “no-w ” vortidty wave 
W= 4 Vorticity wave 2 (downstream only) “no-u” vorticity wave 

It should be noted that while the above more limited wave set is used by the present version of TFaNS, 
CUP3D can handle a more general set of modes by adjusting the number of wave types, so long as 
CUP3D uses the same set of wave types for all acoustic elements. 

We will now discuss how the above more limited wave set is represented in the present version of 
TFaNS. The concepts presented below also apply to CUP3D with a more general set of modes. 

3.2.1 State Vector Conventions 

The acoustic pressure modes for coupling are of the general form: 




Where: 

ipmn ~ mode eigenfunction 

n = harmonic of blade passing frequency 

m = circumferential mode order = nB - kV 
B = number of rotor blades 

V = number of stator vanes 

fi = radial mode order 

Ynty = axial wavenumber normalized by the duct radius 
x =xin Figure 5 but normalized by the duct radius. 


( 3 ) 
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The full expression for perturbation pressure for a given wave type at any interface plane, P, (such as 
upstream going or downstream going waves) is; 


pfyx, r ,<p,t) = /?» 


1 X 


xty-nBQt 


} 


(4) 


The “no— w” vorticity waves have no radial (”w”) velocity component. The defining velocity 
component for this wave is in the axial direction and takes the form: 


u%{r,<p,x,t) = c=o Y 




n = — cok= — ^fi~0 


(5) 


The “ no _ u ” vorticity waves (vorticity with no axial velocity component) are defined with the velocity 
component in the radial direction: 


w p (r,<p,x,t) = Coo Y X 


Y A P n ^W 4u (r)e{^ +9 *^ x "H 


( 6 ) 




n = k= - ™ pt-0 


For coupling purposes these equations are assumed to be valid in the vicinity of each coupling plane, 
P. In these equations all lengths are normalized by the duct radius, fy. The perturbation pressure mode 
amplitudes, A^ nkfi , are normalized by the far-field static pressure, p^ (W=l or 2). The perturbation 

velocity for the vorticity modes (W= 3 or 4) are normalized by the far-field acoustic speed, c OT . The 
acoustic pressure,/?, has the same dimensions as p^. The velocities, u or w, have the dimensions as 
c . The index n counts harmonics of blade passing frequency, m is the circumferential mode order, 
and/* is the radial mode index. The superscript, P, denotes the interface plane where this equation 
is being applied. The subscript, W, denotes the wave type. 

An example of the waves associated with an interface plane is shown in Figure 6 for an isolated stator 
element. Consistent with Figure 1 where interface planes are numbered from inlet to aft, interface 
plane 2 represents the interface plane between the rotor and stator acoustic elements. The interface 
between the stator element and aft radiation element is plane 3. Note that the elements of the state 
vector, A p Wnkfi , shown in this figure have a superscript which denotes the interface plane, P, and a 

subscript which denotes the wave type, W. The source vector elements, B, are the mode amplitudes 
resulting from a prescribed input (such the response of a stator to an incoming wake). 


Interface Plane 


Source Pressure WaveB 



Scattered Pressure Wave A 
Incoming Pressure Wave A 

Incoming Vorticity Wave A 
Incoming Vorticity WaveA 




Source Pressure Wave 
Source Vorticity Wave 
Source Vorticity Wave 
Incoming Pressure Wave 
Scattered Pressure Wccve 

Scattered Vorticity Wave 
Scattered Vorticity Wave 


Figure 6: Mode Scattering through an Isolated Stator 
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To prevent the amplitudes of decaying waves from becoming too large, the axial origin of Equations 
(4) to (6), varies with interface plane and wave type. For a given acoustic element (e.g. Figure 6), 
the origin for each wave type is where each wave exits an acoustic element. For example, in Figure 6, 
the upstream pressure waves, A\, have their origins at interface plane 2 since the upstream pressure 
wave emerges from the element at interface plane 2. Note that A\, for example, is a compressed 
notation collecting all of the state vector elements, A\.. 

Since Equation (4) is the modal representations of acoustic pressure in an annular duct, a form for 
the mode eigenfunction, rp, must be specified. Reference 11 uses the form below for the mode 
eigenfunctions for an annular duct with a mean axial flow. These mode eigenfunctions were originally 
developed by Tyler and Sofrin (Reference 13) and take the form: 

D nm YmjU ( ) 00 

where C mfl and D mfi are constants determined by boundary conditions at the duct walls. For this case 
is the mode eigenvalue and the mode eigenfunction is rpm^ where these modes are normalized 
so that the maximum value of the mode eigenfunction is always a positive one (i.e. xpmax = +i). The 

sJP 

axial wavenumber non-dimensionalized by r d , y Wnkfl , is therefore found to be: 

y, _ " M x (nBM t ) =F - (l - A/|){x: m/< } 2 

y ^nkM = TTW? (8) 

where upstream pressure waves (W = 1) use the upper sign and downstream pressure waves ( W = 2) 
use the lower sign. If Equations (7) and (8) are modified with a kinematic transformation to include 
solid body swirl (as in Reference 2), then Equation (7) does not change. Only Equation (8) changes 
to become: 


^ - M x (nBM, - mM s ) T J(nBM t - mM s f - (l - iW?)jjc^) 2 

Y ^ nkft 1 - M* 

where M s = Swirl Mach Number at the fan tip. 

In the present version of TFaNS (see Reference 2), the “no— w” vorticity wave radial mode function, 
U 3 , from Equation (5) is formulated as a series of cosine waves: 

<V'> ■ -fW) (10) 

To prevent this mode function from becoming infinite when m = 0 in TFaNS, the radial mode function 
in TFaNS uses the form: 


mUi = mcos 


Likewise the “no— u” vorticity wave radial mode function takes on the form: 


( /at{r-r h ) \ 
\ r “- rh ) 


( 11 ) 


mW, 


(r) = ffi co:: H - "A 

“ r r/r^y r d ~r h ) 


( 12 ) 
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The axial wave numbers (non— dimensionalized by rj) associated these waves are. 

~P _ JP _ nBM, — mM s 
Ylnku ~ ^4 nku ~ M x 


( 13 ) 


which corresponds to pure convection of the waves. 

The mode eigenfunctions presented in Equations (7), (9), (11), and (13) are presently used by the 
SOURCE3D source code (Reference 2) as part of the TFaNS Tone FAn Noise design/prediction 
System. The radiation codes utilize Equations (7) and (8) since flow in the inlet and nozzle is swirl 
free. Kousen (Reference 14) developed modes which include arbitrary shear and swirl flows. In the 
future, these modes, or modes like them, will be utilized in a version of the CUP3D code. 

CUP3D assumes that all pressure mode eigenfunctions, no matter how they are formulated, are 
consistent with Equation (4). CUP3D also assumes that all vorticity modes are consistent with 
Equations (5) and (6). In addition all pressure and vorticity modes are normalized so that the 
maximum value of the mode eigenfunction is always a positive one (i.e. VW = +-0- Thus, CUP3D 
is presently capable of accepting other modes (e.g. Reference 14) assuming that the number of wave 
types, W, is finite and a counting system for the wave types can be consistent for all acoustic property 
files entering CUP3D. 

3 2.2 Scattering Coefficient Conventions 

The explanation of the “Kinematics of Modal Scattering” in Reference 1 discusses how modes may 
be scattered into other wave types. Briefly in the CUP3D code: 

Inlets and Nozzles: n, k, fi — » n, k, ft’ 

Stators n, k, [t n, k’, ft’ 

Rotors n, k, fi — » n\ k, fi’ 

The arrow means “is scattered into” so that for example for a stator, an (n, k, fi) mode is scattered into 
a (n, k\ fi’) mode. Therefore stators scatter on circumferential order via A: and radial mode order, fi. 
However, they do not scatter on frequency so that for stators n does not scatter into n’. However, 
rotors are seen to scatter on harmonic via n and radial mode order, fi, while not scattering on k. 

Scattering Coefficients, S, represent the transmission and reflection of waves by an acoustic element. 
An example of mode scattering is shown in Figure 6 for an isolated stator element. This figure was 
discussed in detail in the previous section. Waves entering the stator acoustic element in this example 
are scattered into other waves by the element and then output. For example, the A\ wave can be 
scattered into the A\ wave per the notation in Reference 1 by: 

A 2 «_ c22 A 2 (14) 

Vink' il' - jikii InJqi 

This is the expression for one of the modes (n, k, fi) being scattered by this acoustic element into the 
mode (n, k\ fi’). All other modes are scattered as well within the rules for wave scattering given above 
earlier in this section. 

Note that when mode scattering occurs on harmonic (for a rotor element), noise is scattered into both 
positive and negative harmonics (i.e. — oo < n’ < oo). Thus a positive harmonic, n, may scatter into 
a negative harmonic, n’. This must be taken into account in the scattering matrix output. 

For the purpose of completeness, the total set of equations for a coupled system of inlet radiation, 
rotor, stator, and aft radiation is given by the matrix equation below. This is the same as Figure 4 in 
Reference 1 except for an additional vorticity wave type. Dots indicate blocks of zeros. 
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(15) 


Each S value includes all harmonics (both positive and negative), and circumferential and radial mode 
orders. 

As discussed earlier, this equation can also be represented in a condensed form by Equation 15 in 
Reference 1: 


A = SA +B 


(16) 


Where: 

A = State Vector Array 
S = Scattering Matrix 
B = Source Vector Array 

Equations (4) to (6) suggest that all modes and harmonics are being considered in the calculation. 
However, solution of Equation (15), requires the truncation of the series for a finite size matrix. This 
is partly justified by the fact that acoustic modes are “cut off” above a certain mode order. For the 
acoustic modes given by Equation (7) and Equation (8) cutoff occurs when the the term under the 
square root sign in Equation (8) becomes negative, i.e. 

(nBM,) 2 - (l - Ml){ Kmfl ) \ 2 < 0 (17) 
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Another way of expressing “cut off” is via the cutoff ratio, £, which indicates cut off when £ < 1. For 
modes in a constant area annular duct with a mean axial flow: 


1 = 


nBM t 

Kmfi 


(18) 


Likewise for Equation (9): 


£ = 


nBMj^ - mM s 
Kmfj. 71 — Mx 


(19) 


For the purposes of CUP3D, all the cuton modes up to some sound harmonic, n, should be included 
in the Acoustic Properties File. In addition, a certain number of cutoff modes should be included. 
The issue of how many modes are actually required, however, will be the subject of further study. 
Therefore Equation (4) becomes: 


pfyx,r,<p,t) -p c 


NHT KMAX MUMAX Ix-x^-nBQt] 

111 ^ *4 1 ( 20 ) 


n = -NHT k =KMIN /x = 0 


Where: 


NHT = number of harmonics being run. 

KMIN = minimum value of k being run (different for each harmonic) 

KMAX = maximum value ofk being run (different for each harmonic) 

MUMAX = maximum value offi (different for each harmonic and k-value). 

For vorticity waves Equation (13) indicates that there is no cutoff. This creates a question of how many 
vorticity waves are required for the solution. In the present system, the, (n,k,fi ) , modes are used for 
pressure and vorticity waves. However, this issue will need to be investigated in the future. Also in 
the present SOURCE3D code (Reference 2) only the “no-u” vorticity wave are defined. Thus the 
Equation (5) becomes: 


NHT KMAX MUMAX 

u%(r,<p,x,t ) = Coo 2! 2 2 A ^nku U ^ e 

n = -NHT k = KMIN ft= 0 




( 21 ) 


It should be noted, these Fourier series includes both positive and negative harmonics. For cases 
where tones are not being scattered on harmonic (i.e. for stators and inlet and aft radiation), only the 
positive harmonics need to be used. This is because the negative harmonics are the complex 
conjugates of the positive harmonics. However, for rotor scattering, n’ must include both the positive 
and negative harmonics. Thus the scattering coefficients for the various types of elements take on the 

form: 


Inlet Radiation: 

n > 0, n’ > 0: 

aP' c P'P A p 

A Wnkfi ' °WWnkfi'/Jcfi w nkfi 

(22) 

Rotor: 

n > 0, -NHT< n’ < NHT: 

aP' gP'P A p 

A Wn’kfi’ ^WWn'kti’Nv Wnkfi 

(23) 

Stator: 

n > 0, n’ > 0: 

aP' CP'P A P 

^Wnk'n' ^WWnk'fi'ytkfi tVnkfi 

(24) 

Aft Radiation: 

n > 0, n’ > 0: 

aP' «_ CP'P A P 

A W nkfi' ^W'WnkM'Acpc Wnkfi 

(25) 
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These expressions should be used in formulating the scattering coefficients for CUP3D. Scattering 
coefficients are found by inputting a unit value for A^^ and determining how much is reflected and 
transmitted through the element including all pressure waves and vorticity waves. Thus by setting 
A p Wnktt - 1.0 + iO.O and finding A p Wn , w ^ „ we can determine by using expressions (22) 

through (25). 

The source vector elements, B p Wnktl , are the mode amplitudes resulting from a prescribed outside 
influence such as a rotor wake interacting with a stator vane. B is defined by equations (20) and (21) 
with A w nkfi replaced by B p Wnkf . 
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3.3 REMOVING NEGATIVE HARMONICS FROM CUP3D 

The purpose of this section is to reduce the size of the scattering matrix, S, in Equations (15) and (16), 
by eliminating the need to use negative harmonics in the calculation. The end of the previous section 
pointed out the need to include both positive and negative harmonics in Equation (20). However, for 
fan tones, the negative harmonics are complex conjugates of the positive harmonics. This is the basis 
for the removal of negative harmonics from the calculation. One exception to this as pointed out in 
the previous section, is that the rotor can scatter modes with positive harmonics into modes with 

negative harmonics. This is shown in Equation (23) where can have the n th harmomc 

scatter into the n’th negative harmonic. This aspect of the problem will be included in the discussion 
below. 

We begin with Equation (15) and rewrite it so as to separate the positive harmonics (denoted by a + 
superscript) from the negative harmonics (denoted by a superscript) so that: 



Note that in this equation a scattering coefficient such as S + denotes a positive harmonic, n, 
scattering into a negative harmonic, n’, as is found in a rotor. 

It can be shown that first line in Equation (26) is the complex conjugate of the second line in Equation 
(26). Consequently only one of these lines is needed to calculate a solution. Thus we can use the 

second line: 


A + = S + ~A~ + S + + A + + B + (27) 

Note that all the terms in this equation include only positive harmonic except the term S + ~A~. Also 
note that: 

5 +_ = (S _+ ] a°d A~ = {/4 + } (28) 

so that 

A + = {S" + }*{/4 + }* + S + + /4'*' +B + (29) 

Expanding Equation (29) into real and imaginary parts and separating the real parts and imaginary 
parts of the equation, we get: 



The R subscript denotes the real part of the variable whereas the subscript, 7, denotes the imaginary 
part. Note that in this equation only positive harmonics of BPF are needed as input waves. Also only 
a rotor will create S ~ + coefficients because only a rotor scatters on harmonic. CUP3D uses Equation 
(30) in its calculation of the state vector,/!. 
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3.4 CALCULATION OF POWER LEVELS & FAR-FIELD DIRECTIVITIES 



Figure 7: Calculation of Far— Field Noise From State Vector 

Power levels and far- field directivities are computed from state vector elements at interface planes 
just upstream and just downstream of the noise source. To illustrate the application of this 
computation, the geometry of Figure 1 is utilized. However, CUP3D is not limited to this geometry. 
For the geometry of Figure 1, plane 1 is the upstream interface plane just in front of the fan and plane 
3 is the downstream interface plane just aft of the FEGV. Figure 7 shows only the items directly 
related to the power level and far-field calculations. The state vector elements in Figure 7 are 
denoted by A \ for the upstream pressure waves and A\ for the downstream pressure waves. 

Sound power can be calculated in concept using the equation form below: 

Inlet Power Level, II\ ~ |A}| 2 [P 7 ] Aft Power Level, n\ ~ |A|| 2 [P a ] (31) 

where P[ and P A denote inlet and aft power level multiplication factors, respectively. 

Likewise, far-field pressure directivities can be calculated in concept by applying the equation form 
below: 


Inlet Pressure Directivity, P n ~ [a }][£>,] Aft Pressure Directivity, P n ~ [a|][Z) a ] (32) 

where £>/ and D A are the inlet and aft far-field directivity shapes (based on unit input), respectively. 

The sections which follow will discuss the detailed sound power level and far-field directivity 
equations. Derivations for sound power may be found in Appendix II and derivations for far-field 
pressure directivities may be found in Appendix IV. 
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3.4. 1 Power Level Calculation 

For purposes of the power level calculation only, the following assumptions are made in the vicinity 
of the interface planes where the calculation is taking place. 

• Locally constant area annular duct, 

• Locally uniform axial flow. 

As a result, for the sound power calculation, the standard Tyler/Sofrin modes (Reference 13) are being 
used. A more complicated formulation may be developed in the future. Appendices II and III go into 
further detail on the equations below. 

From Equation (31) and Appendix II, upstream and downstream sound power take on the form. 

fiwnkn ~ 1 A W„kf,\ P Wnkfi 


where n Wn]qi is the non-dimensional sound power for each acoustic mode and harmonic of blade 
passing frequency (BPF). This form of the power is non-dimensionalized by p«,c*.{2r d ) or 
p x Jy a R a Tn (2r d ) 2 . A p Wnkfi is a state vector element where, for the example in Figure 7, P = 1, W = 1 
for upstream pressure waves, and P = 3, W = 2 for downstream pressure waves. The power level 
multiplication factor, Pwniqi ’ is § iven b y : 



nP x I P p 

YaPp V To C 


(1 



+ Mi tp + 


M x 


Rw, 


nkfi\ 


f 

( 


'mu 


(34) 


where the axial Mach Number at plane P, M Xp , is positive downstream and: 

~ P Y Wnkfi 

R Wnkp ~ l \ 

\nBM tp - M x ,y Wnkfl J 


(35) 


Twnkp * s t ^ ie ax i a ^ wavenumber as given by Equation (8) but repeated here as: 

^ ~ M x /iBM !p =F JJnBM,^ 2 - (l - Mj p ) 

yWnkp = 1 - mI 

Xp 

C mu is the radial integral which includes modes of the form given by Equation (7) such that: 

TTlfA. 

1 

| Vliprdr = C md 


(36) 


(37) 
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where r is the radius non-dimensionalized by r^. Using the assumptions above, C mfJL can be integrated 
analytically (see Equation (III. 13) in Appendix III) as: 

( 1 ~ ~ {° 2 ~ (38) 

where rpm^iKm^) and are the tip and hub mode eigenfunction values respectively. Note that 

each mode eigenfunction, xp, is normalized such that xp max = +1- 

The total harmonic sound power level is the sum of the modal powers: 

n'wn- I IKjH* ( 39 > 

k=-°° m=0 

To redimensionalize this non-dimensional power multiply by: p x Jy a RaT (2 r d ) . The units in the 
CUP3D code are: (lbf/in 2 )(ft/sec)(in 2 ) = ft-lbj/sec. To convert to Watts, multiply by 1.3558. Then 
to calculate sound power level, reference this power to 10“ 12 Watts. This is shown as: 



Sound Power = 10 log 


n w ,(p. v5^W 2 ^)(i. 355 8i ^) 


10 - 12 Watts 


(40) 


Where Ra = (53.35 ft lb f /°R lb m )(32.174 lb m ft/sec 2 lb f ) = 1716.5 ft^sec 2 °R and y a = 1.4. Note that 
p 0 o is in lbf/in 2 , is in °R and is in inches. 

3.4.2 Far-Field Directivity Calculation 

The purpose of this section is to present the equations used by CUP3D to calculate far— field 
directivities using the state vector solution (i.e. mode amplitudes) from CUP3D. A derivation of these 
equations may be found in Appendix IV. Appendix V documents the changes which were needed to 
the Eversman inlet and aft radiation codes to make them compatible with CUP3D. Appendix V 
provides an example of how the radiation equations in this section are applied to actual radiation 
codes so they can interface with CUP3D. 

The CUP3D code has been formulated to permit multiple upstream or downstream pressure wave 
types to exist in the duct. This is done in order to accommodate mode sets such as those in Reference 
14. Modes in Reference 14 are formulated in radially dependent axial and swirl flows and are 
therefore more general than the modes we are currently using in TFaNS. The equations below for 
far— field directivities apply to these more general modes thus allowing them to be included in the 
far— field solution. 
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To calculate the far-field directivity two parameters are required: 

1. The solution to the state vector, A p Wnkfi , at the interface plane where a radiation acoustic element 
is connected to the noise source. 

2. Directivity shapes for each propagating pressure wave type and mode from an inlet radiation 
calculation, D Iiy ,(R,d) and an aft radiation calculation, D AWnkfl (R,6). These are determine 
by inputting unit waves into a radiation code and then having the radiation code calculate 
far-field directivities. 


The following far-field pressure directivities are calculated by CUP3D: 

• Total far-field sound pressure level (SPL) directivity for each harmonic 

• Inlet SPL directivity for each harmonic 

• Aft SPL directivity for each harmonic 

• Inlet, aft and total SPL directivities for each circumferential mode order at each 
harmonic (summed over all radial mode orders) 

• Inlet, aft and total SPL directivities for each circumferential and radial mode 
order at each harmonic 

To calculate far-field directivities, start with Equation (IV.24) from Appendix IV, which calculates the 
total far-field pressure directivity for each BPF harmonic: 


Pn(R,0) 





Ankfi 


H = 1 


(41) 


This equation gives the root mean square of the total harmonic pressure pattern averaged 
circumferentially. Circumferential averaging is required so that the azimuthal location of noise 
measurement relative to a circumferential origin on the engme is not needed to perform the 
calculation. Multiple circumferential mode orders will cause azimuthal variations in the far field 
directivity which are averaged by this equation. 


To calculate the far-field SPLs, note thatP„(i?,0) are normalized by the far-field static pressure, p^. 
Therefore: 


where P re f = 20pPa. 


SPL nml = 20.og(^^) 


(42) 



In Equation (41), p lnkfi (R,0 ) refers to the inlet complex pressure directivity for each (m,p) mode at 

the interface plane, P, where an inlet is connected to the noise source (e.g. in Figure 7, P = 2). This 
is calculated as: 


PlnktiW’Q) 


^UfSTR 

X, A WnlqPlWnkfi ^ ) 


W= 1 


(43) 


Where: 

A p Wnkfl = Upstream propagating pressure wave amplitudes from the state vector at interface plane, P 
D IWnkfi {P, 0) = Mode inlet directivity shape functions 

CUP3D is configured to use multiple upstream pressure wave types so that these pressures are added 
together to obtain the complex pressure directivity for each (m, p) mode, p ]nlqi (R,6). The number 

of forward propagating pressure wave types is given as N UPS tr. Note that for all modal representations 
developed to date (e.g. Reference 14) only one type of forward propagating mode has been found (i.e. 
Nupstr = !)• However, CUP3D uses Equation (43) as shown to permit more upstream pressure wave 
types to be added, if necessary. 

In Equation (41), refers to the aft complex pressure directivity for each (m,p) mode at 

the interface plane, P, where the nozzle is connected to the noise source (e.g. in Figure 7, P = 3). This 
is calculated as: 


n presdir 

PAnkfi(R>0) = X A WnkfPAWnkfi( R ’ 0) (44) 

w=n ups ra + 1 

Where: 

A p WnJqjl = Downstream propagating pressure wave amplitudes from the state vector at interface plane, P 
D AWnktl {R,6) = Mode aft directivity shape functions 

CUP3D is configured to use multiple downstream pressure wave types so that these pressures are 
added together to obtain the complex pressure directivity for each (m,p) mode, p Ankfl (R,0). The 

number of propagating pressure wave types is assumed below to be Npresdir ■ For the modes presently 
in TFaNS and discussed in Section 3.2, Npresdir = 2. 

Using Equation (41), the other types of directivities calculated by CUP3D can be formulated. To 
calculate the total inlet directivity (circumferentially averaged), set p Ankfl = 0 in Equation (41). T his 

yields the summation over all inlet modes to get the inlet pressure directivity where: 


V *=-00 


X &p 


Inkfi 


lu=0 


(45) 
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The same can be done with aft directivity where: 


JV<K.»>= / X 

Y -00 

To determine total, inlet, and aft pressure directivities for each circumferential mode. Equations (41), 
(45), and (46), respectively, may be applied without the summation over the circumferential mode 
orders to get the root mean square pressure directivity for each circumferential mode order. Thus 
total, inlet, and aft circumferential mode pressure directivities become: 


Ztip. 


Ank/u 


\fi=0 


(46) 




X ftp ** + X ft*- 


Ankfi 


Lu=0 


/i = 0 


P+uJM) 



PnkaftW’V 



(47) 


(48) 


(49) 


Note that these results are not circumferentially averaged since there is no circumferential variation 
of a single circumferential mode order. 

The same process may be used to get total, inlet, and aft pressure directivities for particular (m,u) 
modes Where equations (47), (48), and (49) become respectively: 


P n ~ y | fiPlnkn + ^PAnkfi\ 

(50) 

P nkfi in i et ^ R ^ = VWl 

(51) 

R nkfi a jf R ’^ = -J\ftPAnkfi\ 

(52) 


SPEs can then be calculated for these cases by inserting these pressure directivities into Equation (42). 
CUP3D uses Equations (41) to (52) to calculate the various far-field directivities. 
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4. CONCLUDING REMARKS 

A Tone Fan Noise Design/Prediction System (TFaNS) was developed. The purpose of this system is 
to predict tone noise emanating from a fan stage including the effects of reflection and transmission 
by the rotor and stator and by the duct inlet and nozzle. These effects have been added to an existing 
annular duct/isolated stator noise prediction capability. 

TFaNS consists of: 

• The codes that compute the acoustic properties (reflection and transmission coefficients) of the 
various “acoustic elements” (e.g. inlet, rotor, stator, nozzle). These properties are written to 
Acoustic Properties Files, 

• CUP3D: Fan Noise Coupling Code that reads these files, solves the coupling problem, and outputs 
the desired noise predictions, 

• AWAKEN: CFD/Measured Wake Postprocessor which reformats CFD wake predictions and/or 
measured wake data so they can be used by the system. 

This report discussed the organization of TFaNS and then provided the technical documentation for 
the CUP3D Fan Noise Coupling Code. It also documented for code developers the file format for 
Acoustic Properties Files required as input to CUP3D. The Eversman radiation codes (which are 
presently part of TFaNS) were made available as an example to show how information for the Acoustic 
Properties Files could be created. This last portion of the report also acted to document the changes 
to the Eversman codes in order to make them part of TFaNS. 
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LIST OF SYMBOLS 


PARAMETERS 

A — State Vector 

Awnkn * $ tate V ector Element 
B = Source Vector 

Bwnkfi = S° urce Vector Element 

B = number of blades 

c — acoustic speed 

D IiVnkfi = inlet directivity from an inlet radiation code at radius, R in the far— field given unit input. 

D Awnkn = a ft directivity from an aft radiation code at radius, R in the far— field given unit input. 

I = instantaneous acoustic intensity 

I = time average intensity of I 

J — Bessel Function of the first kind 

k = circumferential mode number from m — nB — kV where — oo < k < +oo 
kx = axial wavenumber in the Eversman radiation codes 

k x — axial wavenumber non —dimensionalized by r^ in the Eversman radiation codes 

m — circumferential mode order where m is positive in the direction of rotor rotation 

Mf — mean axial Mach Number at the input plane ( interface plane) of the Eversman radiation codes 

M s — tip swirl Mach number 

M t = tip rotational Mach number 

M x = mean axial Mach number 

n = harmonic of blade passing frequency 

Nj = uncorrected fan speed (rpm) 

N] c = corrected fan speed (rpm), see Equation (V5) 

Nh = number of BP F harmonics being calculated (see Appendix II) 

N p = number of interface planes 

p = acoustic pressure 

Poo = reference pressure (static far— field pressure) 
p P = mean static pressure at interface plane, P 

Pwnkp d ~ m °d a l acoustic pressure in the Eversman radiation codes 

Pnkp (R, 8) = far— field pressure directivity for a given mode and harmonic non —dimensionalized by p^ 
Pinkfi = inlet f ar ~ field pressure directivity for a given mode and harmonic non— dimensionalized by p^ 

PAnkfi = a ft far— field pressure directivity for a given mode and harmonic non —dimensionalized by p^ 

P(R, 8, <p, t) — time dependent pressure at a given point in the far— field non —dimensionalized by p^ 
P n (R,8) = Harmonic far— field pressure directivity non— dimensionalized by p^ 
r = radial direction non —dimensionalized by r d 

r = dimensional radius (inches) 

r<i = Outer Duct Radius (inches) 

rh = Hub Duct Radius (inches) 

R = Radial distance from the engine origin to some location in the far— field. 

Ra =Idealgas constant for air = (53.35 ft lbf/°R lb m )(32.174 lb m ft/sec 2 Ibf) = 1716.5 ft 2 /sec 2 °R 
Rwnkp = relates acoustic velocity modes to acoustic pressure modes as in equation (II. 5) 

RWnfy = relates acoustic velocity modes to acoustic pressure modes (non-dimensional), see eqn. (II.6) 
S = Scattering matrix 
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= Scattering Coefficient where “primed’’ subscripts denote the wave which is scattered into, 
time (sec) 

mean static temperature (°R) 

total temperature in the far field or in the inlet (°R) 
acoustic velocity 

acoustic velocity in the axial direction 

mean flow axial velocity . 

Radial Mode Shape (mode eigenfunction) for the axial component of the acoustic velocity 

number of vanes 
swirl flow velocity 

acoustic velocity in the radial direction 

Radial Mode Shape (mode eigenfunction) for the radial component of the acoustic velocity 
axial direction (positive downstream) non-dimensionalized by r^ 

- axial direction (dimensional) 

- Bessel Function of the second kind 

- axial wavenumber 

= axial wavenumber non-dimensionalized by rj 
= ratio of specific heats = 1.4 for air 
= constant used to calculate solid body swirl 

= Reduced frequency in the Eversman Radiation codes, see Equation (VI) 

= directivity angle of the far -field location where data is to be measured. 

= [ n B - kVyp - nBQt (See Appendix IV) 

= unknown constant in Bessel’s equation (Appendix III) 

= Mode eigenvalue non-dimensionalized by rj 
= radial mode order 
= acoustic sound power level 

= acoustic sound power level non—dimensionalized by p 00 Coo(2r (i ) 

= mean flow static density 
= = hub to tip ratio of the duct 

= circumferential direction (positive in the direction of rotor rotation) 

= perturbation velocity potential mode amplitude in the Eversman radiation codes 


Wwnkn = P ertur b at i° n velocity potential for a single mode in the Eversman radiation codes 
Vwnkfi^’ ^ = f ar ~fi e ^ perturbation velocity potential mode directivity in the Eversman radiation codes 

rp = Radial Mode Shape (mode eigenfunction) for the pressure waves normalized by xpmax 
Vmax = Maximum value of each radial mode shape (mode eigenfunction) which is defined as +1 
a) = rotational frequency (rad/sec) 

Q = rotor rotational speed (rad! sec) 


25 



SUBSCRIPTS AND SUPERSCRIPTS 

aft = Aft radiated direction 

f = At the input plane (interface plane) of the Eversman radiation codes 

inlet = Inlet radiated direction 

k = circumferential mode number from m = nB - kV where — oo < k < +oo 

m = circumferential mode order where m is positive in the direction of rotor rotation 

n — harmonic number 

P = Interface Plane Number 

W = Wave type, e.g. for modes in a uniform axial flow with solid body swirl: 

= 1 (upstream pressure), 2 (downstream pressure), 3 &4 (downstream vorticity) 
p = radial mode order 

oo = reference value (i.e. in the far— field) 
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APPENDIX I: CREATING A CUP3D ACOUSTIC PROPERTIES 

FILE 

1.1 ACOUSTIC PROPERTIES FILE GENERAL INFORMATION 

Acoustic Properties Files are output by rotor, stator, inlet radiation and aft radiation computer codes 
for input into CUP3D. An Acoustic Properties File is required for each reflecting acoustic element 
in the system. No Acoustic Properties File is needed for a non— reflecting element. 

For a code to create an Acoustic Properties File, the following information should be placed in a file 
in the order shown below. Note that the code generating this file is designated by a specific type 
(ICODE). Based on the type of code, the information required by CUP3D will vary. 

The parts of the Acoustic Properties File are: 

• Header Cards: Contain information about the run including the code which created the file, the 
acoustic modes contained in the file, and certain geometry and performance of the acoustic 
element. This information is split into information which is common to all elements and 
information which is common to a specific type of element. 

• Scattering Coefficients: These are the reflection and transmission coefficients of the acoustic 
element. Definitions for these coefficients are given in Sections 3.2 along with a method for 
creating them. 

• Source Vector The source vector is output if the acoustic element is influenced by an outside force 
such as a wake interacting with a vane or blade. In this case acoustic mode amplitudes and vorticity 
waves are created which form the source vector and are placed in the Acoustic Properties File. 
Definitions for these coefficients are given in Sections 3.2. 

• Far— field Directivities: If the acoustic element radiates noise to the far— field, then far-field 
directivities are placed in the Acoustic Properties File. Pressure directivities result from unit input 
for each duct mode into a radiation code. These directivities are given at a constant radius, R, (See 
Figure 5) from a user specified origin. 

Multiple cases can be placed in the same file. In this situation, the first case contains all of the cards 
given below as defined by the type of acoustic element being considered. All of the cases which follow 
will then omit Card 1. The cases are then stacked one below the other. See stator example 2 in Section 

1.3.7. 

The following terminology is utilized by the system (see Figure 3 and Figure 4): 

• An acoustic element (or region) refers to a part of the system with acoustic properties such as a 
rotor, stator, inlet radiation or aft radiation. An Acoustic Properties File must be created for any 
acoustic element with acoustic reflection properties. An Acoustic Properties File is not needed 
for a non -reflecting region. 

• An interface plane is the plane where two acoustic elements meet. 

• Scattering coefficients correspond to reflection and transmission coefficients for a particular 
region (see Section 3.2). 

• The Source vector corresponds to the mode amplitudes resulting from a prescribed outside 
influence such as a rotor wake interacting with a stator vane. 
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Variables in this file are non— dimensionalized using the following conventions: 

• Lengths: Non -dimensionalized by fan tip diameter at fan leading edge (inches), DFAN. 

• Temperature: Non— dimensionalized by the far- field static temperature (°R) (far from the 
nacelle), TINF. 

• Pressure: Non -dimensionalized by the far-field static pressure (psia) (far from the nacelle), 
PINF. 

• Vorticity : Non— dimensionalized by c* where c x = Jy a R a {TINF} 

1.2 ACOUSTIC PROPERTIES FILE ASSUMPTIONS 

• Waves are assumed to be given to the code in the following order: upstream propagating waves 
first, then downstream propagating waves. Pressure waves are input before vorticity waves. 
Vorticity waves only propagate in the downstream direction. NWAVETYPE denotes the total 
number of waves which can propagate. NUPSTRE = number of upstream propagating pressure 
waves. NPRESDIR = the total number of propagating pressure waves. In the present code: 

NWAVETYPE = 4, NUPSTRE = 1, NPRESDIR = 2. 

• All source vector information in the Acoustic Properties Files should be expressed as peak 
complex modal amplitudes. 

• The code assumes that all acoustic elements were created with the same number of blades 
(NBLADE), vanes (NVANE), and fan rotational tip speed (VTIPC). Hub to tip ratios at the 
interface planes should be the same for the two connecting acoustic elements for the modes to be 
passed correctly between regions. The code, however, does not prevent hub to tip ratios from 
differing from region to region. 

• The code presently assumes that all acoustic elements are in a single duct (such as in Figure 1). 
The concept of a core engine duct and bypass duct is not implemented in the code. 

• The code is written so as to assume that all cases are being run in a user defined order by CUP3D. 
The same order is assumed for all acoustic elements. Therefore, if multiple cases were run by a 
stator code for example using a different order than was used for a rotor code, then CUP3D will 
not run correctly. 

• Far field angles are assumed to be stored in the Acoustic Properties Files with the smallest angle 
first and the largest angle last. Zero degrees is assumed to be in front of the inlet. 

• It is assumed that the user has chosen an origin for the entire calculation and that the origin is 
consistent for all Acoustic Properties Files being input into CUP3D. Far-field directivities 
should utilize this origin to specify distance from the engine and angles around the engine. 

1.3 ACOUSTIC PROPERTIES FILE FORMAT 

The discussion of the Acoustic Properties File format will be as follows. First the header cards 
common to all cases will be discussed in Section 1.3.1. After that all other cards will be discussed based 
on the type of computer code which is outputting this file (i.e. inlet radiation code in Section 1.3.2, 
rotor code in Section 1.3.3, stator code in Section 1.3.4, aft radiation code in Section 1.3.5). Details 
of the scattering coefficient and source vector formats is found in Section 1.3.6. Examples of Acoustic 
Properties Files are found in Section 1.3.7. Also the method used to modify the Eversman radiation 
codes is given in Section V.2, Appendix V. Appendix V may be helpful towards the understanding of 
the process for creating files. 
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The input data file structure has distinct blocks of data referred to as cards. Each card begins on a 
new line and the input data is either free format or formatted as given in the card. 

Note that a ^ denotes parameters which must have the same value from element to element for the 
CUP3D code to run correctly. It is imperative that these parameters be defined exactly as given in 

this document. 

Except where noted, all input is free format. 

For all complex numbers the following format statement is used: 

999 FORMAT(10E13.6) 

1.3. 1 Header Cards for all Acoustic Properties Files 

All Acoustic Properties Files contain the header cards below. 

Card 1: (for multiple cases this card is output only once with the first case) 

• ICODE: Computer Code which created this file 

Note: ICODE is designed so that the first number denotes the type of code which created this 
file. The second number refers to which code created this file. For example: ICODE = 11: 

The first 1 refers to inlet radiation codes, the second 1 refers to the Eversman inlet radiation 

code. 

= 10 Other Inlet Radiation Code 
=11 Eversman Inlet Radiation Code 
= 12 Caruthers Inlet Radiation Code 
=20 Other Aft Radiation Code 
=21 Eversman Aft Radiation Code 
=22 Caruthers Aft Radiation Code 
=30 Other Rotor Code 
=31 SOURCE3D Rotor Code 
=40 Other Stator Code 
=41 SOURCE3D Stator Code 

• NCASE: Number of “cases” where a case is defined as all harmonics for a given condition or 
configuration. Acoustic Properties Files are read in the order ICASE = 1, NCASE (see Card 2). 

Card 2: ICASE: Case Number: FORMAT(lX,’Case Number’, 13) 

Card 3: CODETITLE: Title Card - identifying the case. Use the title from the code which was run. 
FORMAT(A80). 
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Card 4: 


• NHT: Number of blade passing frequency harmonics run. 
ifc VTTPC: Corrected fan rotational tip speed (ft/sec), 
where: VTIPC = ^N lc ( DFAN ) 

DFAN = Fan tip diameter (inches) at the fan leading edge 


N lc = Corrected Fan Rotational Speed (rpm) 



Nj = Fan rotational speed (rpm) 

T 0 i = reference temperature (°R) used to correct the fan speed to a “standard day”. 
For zero Mach number in the far -field, this temperature is the static 
far— field ambient temperature, TINF. 


Card 5: 


^ NBLADE: Number of rotor blades 
^ NVANE: Number of stator vanes. 

Card 6: (KMIN(N,IWAVE), KMAX(N,IWAVE), N=1,NHT), IWAVE = 1 ,NWAVETYPE 
Where for the case which was run: 

KMIN = Minimum circumferential mode number, K, run for each harmonic, 

KMAX = Maximum circumferential mode number, K, run for each harmonic where: 
m = n* NBLADE - K* NVANE 
n = harmonic number 
m = circumferential mode order. 

Card 7: (((M(N,K,IWAVE), K=KMIN(N, IWAVE), KMAX(N, IWAVE) ), N=1,NHT), IWAVE = 
1,NWAVETYPE): Circumferential mode order, m, where m is positive in the direction of rotor 
rotation. 


Card 8: (((MUMAX(N,K,IWAVE), K=KMIN(N, IWAVE), KMAX(N, IWAVE) ), N=1,NHT), 
IWAVE = 1,NWAVETYPE): Largest Value for fi for each circumferential mode order and 
harmonic, fi refers to the radial mode order where 0 </u < MUMAX. 

Card 9: WRITE(IUNIT,999) ((((GAMMA(N,K,MU,IDIR), MU=0,MUMAX(N,K,IDIR)), 
K=KMIN(N,IDIR), KMAX(N,IDIR)), N=1,NHT), IDIR = 1,NPRESDIR): Complex Axial 
wave numbers for acoustic pressures directed upstream (IDIR=1) or downstream (IDIR=2) 
normalized by the duct radius. This parameter is presently input as one set of values per region. 
It is anticipated that at some time in the future, separate upstream and downstream GAMMA’s 
will be input for each interface plane. 
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1.3.2 Inlet Radiation Code: 10 < ICODE < 20 
I.3.2.1 Header Cards 

Information in Cards 11 through 16 are only used to confirm for the user which Acoustic Properties 
File the user has chosen. These parameters are not used in code calculations. While parameters in 
Card 10 are not presently required for code calculations, this information expected to be required for 
code calculations in the future. 

Card 10: 

• SIGMAR: Hub to tip ratio of the duct at XINTER. This number should, strictly speaking, be the 
same for two codes meeting at a given interface plane. 

• XINTER: Axial location of the interface plane relative to the origin (non-dimensionalized by 
DFAN) as specified by the user. 

• MINF: Far-field axial Mach number (far from the nacelle, positive downstream) 

• MDUCT: Mass averaged axial Mach number at the interface plane (positive downstream) 

Card 11: NLINO: Number of acoustically lined elements (outer wall). Set this value equal to zero if 
there is no outer wall liner. 

If ( NLINO .GT. 0 ) THEN 

Card 12: 

• MBEGO: Outer wall element number for start of liner. 

• XBEGO: Axial location where outer liner begins relative to the origin (non-dimensionalized by 
DFAN) as specified by the user. 

• XENDO: Axial location where outer liner ends relative to the origin (non-dimensionalized by 
DFAN) as specified by the user. 

Card 13: WRITE(IUNIT,999) ZIMPO(I),I=MBEGO,NLINO+MBEGO: Liner impedance 

(Resistance, Reactance) for each element on outer wall non— dimensionalized by local pc (i.e. 
density x acoustic speed). 

ENDIF 

Card 14: NLINI: Number of acoustically lined elements (inner wall). Set this value equal to zero if 
there is no inner wall liner. 

If ( NLINI .GT.0) THEN 

Card 15: 

• MBEGI: Inner wall element number for start of liner. 

• XBEGI: Axial location where inner liner begins relative to the origin (non-dimensionalized by 
DFAN) as specified by the user. 

• XEN DI: Axial location where inner liner ends relative to the origin (non-dimensionalized by 
DFAN) as specified by the user. 
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Card 16: WRITE(IUNIT,999) ZIMPI(I),I = MBEGI,NLINI + MBEGI: Liner impedance 

(Resistance, Reactance) for each element on the inner wall non— dimensionalized by the local 
gc (i.e. density x acoustic speed). 

ENDIF 

1.3. 2. 2 Scattering Coefficients (see Section 1.3.6) 

Card 17: Write: FORM AT(1X, ’Scattering Matrix and Source Vector’) 

Card 18: Write: FORMAT(lX, ’Scattering Matrix’) 

Card 19: Scattering coefficients are written as follows: 

C 

C WRITE SCATTERING MATRIX INFORMATION, S 21 , AT THE INTERFACE 
C PLANE DOWNSTREAM OF THIS REGION. THESE SCATTERING 
C COEFFICIENTS CORRESPOND TO UPSTREAM PROPAGATING WAVES 
C SCATTERED INTO DOWNSTREAM PROPAGATING WAVES. NOTE THAT 
C S 31 & S 41 ARE NOT WRITTEN SINCE THEY ARE VORTICITY WAVE BASED. 

C THERE ARE NO VORTICITY WAVES EMANATING FROM THE INLET 
C RADIATION CODE. THUS S 31 AND S 41 ARE ZERO. 

C 

DO 10 IW = 1,NUPSTRE 
DO 20 IWP = NUPSTRE + 1 ,NPRESDIR 
C 

WRITE(IUNIT,999)(((( S^N, K, MU, MUP), 


& 

MUP = 

0,MUMAX(N,K,IWP)), 


& 

MU = 

0,MUMAX(N,K,IW)), 


& 

K = KMIN(N,IW),KMAX(N,IW)), 

<-k 

& 

c 

N =1, 

NHT) 

<— n 

20 

CONTINUE 



10 

CONTINUE 



i.3.2.3 

Source Vector (see 

Section 1.3.6) 


Card 20: IFORCE = 0 No Source Vector is Input 



= 1 Source Vector is input 


IF IFORCE = 1: 

Card 21: Write: FORMAT(lX, ’Source Vector’) 


Card 22: Upstream propagating pressure mode amplitudes normalized by PINF 

DO 10 IW = 1, NUPSTRE 
READ(IUNIT,999) ((( MU,K,N), 

& MU =0, MUMAX(N,K,IW)), 

& K = KMIN(N,rW),KMAX(N,TW)), 

& N = 1, NHT) <- n 

10 CONTINUE 
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Card 23: Downstream propagating pressure mode amplitudes normalized by PINF 


DO 10 IW = NUPSTRE + 1,NPRESDIR 
READ (IUNIT,999) ((( fl£(MU,K,N), 

& MU = 0, MUMAX(N,K,IW)), 

& K = KMIN(N,IW),KMAX(N,IW)), <-* 

& N =1, NHT) 

10 CONTINUE 

Card 24: Downstream propagating vorticity waves normalized by c^. 


DO 10 IW = NPRESDIR + 1, NWAVETYPE 
READ(IUNIT,999) ((( B%( MU,K,N), 

& MU = 0, MUMAX(N ,K,IW)), 

& K = KMIN(N,IW),KMAX(N,IW)), <-* 

& N =1, NHT) 

10 CONTINUE 


ENDIF 

13.2.4 Far-field Directivities 

Card 25: Write: FORMAT(lX,’Farfield Directivities’) 


Card 26: 

• NANGLE: Number of polar angles ( 6 in Figure 5) in the directivity pattern. The code will 
interpolate these angles to form an array every 1° from 0°s ANGLE < 180°. In the present code, 
any angle which is within 25° of the aft most computational boundary will be discarded. This is 
consistent with the baffle boundary condition which presently is found in the Eversman Inlet 
Radiation code. 

• DARRAY: Radial distance ( R in Figure 5) from the origin non-dimensionalized by DEAN 
(origin specified by the user) of the polar directivity array. 

Card 27: Far-field polar angle in degrees (0° is in front of the inlet) 

• WRITE(IUNIT, * ) ( ANGLE(IANGLE), IANGLE = 1, NANGLE) 

Card 28: Write: FORMAT(lX, ’COMPLEX PRESSURE DIRECTIVITIES’) 

Card 29: Complex pressure directivities (real and imaginary parts of the directivities for each mode). 


DO 10 IW = 1, NUPSTRE 

WRITE(IUNIT,999) (((( DIRECTIN (N,K,MU ,IANGLE), 
IANGLE = 1, NANGLE), 

MU = 0, MUMAX(N,K,IW)), 

K = KMIN(N,IW), KMAX(N,IW)) 

N =1, NHT) 

10 CONTINUE 
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1.3.3 Rotor Code: 30 < ICODE < 40 
I.3.3.1 Header Cards 

Information in Cards 12 through 16 are only used to confirm for the user which Acoustic Properties 
File the user has chosen. These parameters are not used in code calculations. Parameters in Cards 
10 and 11 (with the exception of XI and X2) are used to calculate power in the code. XI and X2 are 
not presently used for calculations in the code, but are expected to be needed in the future. 

Card 10: Geometry and performance at the interface plane upstream of the region 

• SIGMAR1: Hub to tip ratio of the duct at XI. This number should, strictly speaking, be the same 
for two codes meeting at a given interface plane. 

• XI: Axial location of the interface plane upstream of the acoustic element relative to the origin 
(non— dimensionalized by DFAN) as specified by the user. 

• DLA1: Outer duct diameter of the interface plane upstream of the acoustic element 
(non— dimensionalized by DFAN). 

• MX1: Mass averaged axial Mach number at the interface plane upstream of the acoustic element. 

• MTIP1: tip rotational Mach number at the interface plane upstream of the acoustic element. 

• PS1: Mass averaged static pressure at the interface plane upstream of the acoustic element 
(non— dimensionalized by PINF). 

• TS1: Mass averaged static temperature at the interface plane upstream of the acoustic element 
(non— dimensionalized by TINF). 

Card 11: Geometry and performance at the interface plane downstream of the region 

• SIGMAR2: Hub to tip ratio of the duct at X2. This number should, strictly speaking, be the same 
for two codes meeting at a given interface plane. 

• X2: Axial location of the interface plane downstream of the acoustic element relative to the origin 
(non— dimensionalized by DFAN) as specified by the user. 

• DIA2: Outer duct diameter of the interface plane downstream of the acoustic element 
(non -dimensionalized by DFAN). 

• MX2: Mass averaged axial Mach number at the interface plane downstream of the acoustic 
element. 

• MTIP2: tip rotational Mach number at the interface plane downstream of the acoustic element. 

• PS2: Mass averaged static pressure at the interface plane downstream of the acoustic element 
(non— dimensionalized by PINF). 

• TS2: Mass averaged static temperature at the interface plane downstream of the acoustic element 
(non— dimensionalized by TINF). 

Card 12: NDAT: Number of streamlines where geometry is specified 

Card 13: DLA(I), 1=1, NDAT: Streamline diameter (non— dimensionalized by DFAN) 

Card 14: RSOLIDITY(I), 1=1, NDAT: Rotor solidity (rotor chord/rotor pitch) 

Card 15: RALPHCH(I), 1=1, NDAT: Rotor stagger angle relative to the circumferential direction 
(degrees). This is normally considered to be a number between 0 and 110 degrees. 
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Card 16: SSOC(I), I=1,NDAT: Upstream stator trailing edge to rotor leading edge aerodynamic 
Spacing to upstream stator aerodynamic chord for each streamline (Set equal to zero if there 
is no upstream stator cascade). 

I.3.3.2 Scattering Coefficients (see Section 1.3.6) 

Card 17: Write: FORMAT(lX, ’Scattering Matrix and Source Vector’) 

Card 18: Write: FORMAT(lX, ’Scattering Matrix’) 


Card 19: Scattering coefficients as follows 


C 

c 

c 

C 

C 

c 

c 


WRITE SCATTERING MATRIX INFORMATION S 12 , S I3 , S 14 , AT THE 
INTERFACE PLANE UPSTREAM OF THIS REGION. THESE SCATTERING 
COEFFICIENTS CORRESPOND TO DOWNSTREAM PROPAGATING WAVES 
ENTERING THE UPSTREAM INTERFACE PLANE AND BEING REFLECTED 
BY THE ELEMENT BACK OUT AS UPSTREAM WAVES. 


DO 30 IW = NUPSTRE+ 1,NWAVETYPE 
DO 40 IWP = 1, NUPSTRE 


WRITE(IUNIT,999) (((( ( NP, K, MU, MUP), 


& 

MU = 0, MUMAX(N,K,IW)), 


& 

MUP = 0, MUMAX(LABS(NP),K,IWP)), 


& 

K = KMIN(N,IW),KMAX(N,IW)), 

k 

& 

N = 1, NHT), 

<r- ft 

& 

NP = -NHT, NHT) 

n 9 

40 

CONTINUE 


30 

CONTINUE 



C 

C 

c 

c 

c 

c 

c 


WRITE SCATTERING MATRIX INFORMATION, Su, AT THE INTERFACE 
PLANE UPSTREAM OF THIS REGION. THESE SCATTERING COEFFICIENTS 
CORRESPOND TO UPSTREAM PROPAGATING WAVES AT THE 
DOWNSTREAM INTERFACE BEING TRANSMITTED THROUGH THE 
ACOUSTIC ELEMENT TO THE UPSTREAM INTERFACE 

DO 50 IW = 1, NUPSTRE 
DO 60 IWP = 1, NUPSTRE 


WRITE(IUNIT,999)((((( Sj^N, NP, K, MU, MUP), 

& MU = 0, MUMAX(N,K,IW)), 

& MUP = 0, MUMAX(LABS(NP),K,IWP )), 

& K = KMIN(N,IW),KMAX(N,IW)), «-* 

& N =1, NHT), 

& NP = -NHT, NHT) n 


C 

60 CONTINUE 
50 CONTINUE 
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c 

C WRITE S 22 > $32> S 42 , S 23 , S 33 , S 43 , S 24 , S 34 , S 44 . THESE COEFFICIENTS 
C CORRESPOND TO DOWNSTREAM PROPAGATING WAVES INPUT INTO THE 

C REGION AT THE UPSTREAM INTERFACE PLANE WHICH ARE 

C TRANSMITTED THROUGH THE ELEMENT TO THE DOWNSTREAM 

C INTERFACE PLANE 
C 

DO 50 IW = NUPSTRE+ 1, NWAVETYPE 
DO 60 IWP = NUPSTRE + 1 , NWAVETYPE 


WRITE(IUNIT,999)((((( S^N, NP, K, MU, MUP), 


& 

MU =0, MUMAX(N,K,IW)), 


& 

MUP = 0, MUMAX(IABS(NP),K,IWP)), 


& 

K = KMIN(N,IW),KMAX(N,IW)), 


& 

N = 1, NHT), 

<■— n 

& 

NP = -NHT, NHT) 

<— n 


60 CONTINUE 
50 CONTINUE 

C 

C WRITE SCATTERING MATRIX INFORMATION, S 2 i, S 31 , S 4h AT THE 
C INTERFACE PLANE DOWNSTREAM OF THIS REGION. THESE SCATTERING 
C COEFFICIENTS CORRESPOND TO UPSTREAM PROPAGATING WAVES 
C REFLECTED BY THE ACOUSTIC ELEMENT AS DOWNSTREAM 
C PROPAGATING WAVES. 

C 


DO 70 IW = 1, NUPSTRE 
DO 80 IWP = NUPSTRE + 1, NWAVETYPE 

WRITE(IUNIT,999)((((( SgP^N, NP, K, MU, MUP), 


& 

MU =0, 

MUMAX(N,K,IW)), 


& 

MUP = 0, 

MUMAX(IABS(NP),K,IWP)), 


& 

K = KMIN(N,rW),KMAX(N,IW)), 

<-Jt 

& 

N =1, 

NHT), 

<— n 

& 

NP = -NHT, 

NHT) 

<— n 

1 

1 

CONTINUE 

CONTINUE 
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I.3.3.3 Source Vector (see Section 1.3.6) 

Card 20: IFORCE = 0 No Source Vector is Input 
= 1 Source Vector is input 


IF IFORCE = 1: 

Card 21: Write: FORMAT(lX, ’Source Vector’) 

Card 22: Upstream propagating pressure mode amplitudes normalized by PINF 

DO 10 IW = 1,NUPSTRE 
READ(IUNIT,999) ((( B J(MU,K,N), 

& MU =0, MUMAX(N,KJW)), 

& K = KMIN(N,IW),KMAX(N,IW)), 

& N = 1, NHT) 

10 CONTINUE 

Card 23: Downstream propagating pressure mode amplitudes normalized by PINF 


<r-k 

<— n 


DO 10 IW = NUPSTRE + 1,NPRESDIR 
READ (IUNIT,999) ((( flg<MU,K,N), 

& MU = 0, MUMAX(N,K,IW)), 

& K = KMIN(N,IW),KMAX(N,IW)), «"* 

& N = 1, NHT) *~ n 

10 CONTINUE 

Card 24: Downstream propagating vorticity waves normalized by c^. 


DO 10 IW = NPRESDIR + 1, NWAVETYPE 
READ(IUNIT,999) ((( R^(MU,K,N), 

& MU = 0, MUMAX(N,K,IW)), 

& K = KMIN(N,IW),KMAX(N,IW)), 

& N =1, NHT) <~ n 

10 CONTINUE 

ENDIF 
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1.3.4 Stator Code: 40 < ICODE < 50 
1.3.4. 1 Header Cards 

Information in Cards 12 through 16 are only used to confirm for the user which Acoustic Properties 
File the user has chosen. These parameters are not used in code calculations. Parameters in Cards 
10 and 11 (with the exception of XI and X2) are used to calculate power in the code. XI and X2 are 
not presently used for calculations in the code, but are expected to be needed in the future. 

Card 10: Geometry and performance at the interface plane upstream of the region 

• SIGMAR1: Hub to tip ratio of the duct at XI. This number should, strictly spe akin g, be the same 
for two codes meeting at a given interface plane. 

• XI: Axial Location of the interface plane upstream of the acoustic element relative to the origin 
(non - dimensionalized by DFAN) as specified by the user. 

• DIAl: Outer duct diameter of the interface plane upstream of the acoustic element 
(non— dimensionalized by DFAN) 

• MX1: Mass averaged axial Mach number at the interface plane upstream of the acoustic element. 

• MTIP1: tip rotational Mach number at the interface plane upstream of the acoustic element. 

• PS1: Mass averaged static pressure at the interface plane upstream of the acoustic element 
(non— dimensionalized by PINF). 

• TS1: Mass averaged static temperature at the interface plane upstream of the acoustic element 
(non— dimensionalized by TINF). 

Card 11: Geometry and performance at the interface plane downstream of the region 

• SIGMAR2: Hub to tip ratio of the Duct at X2. This number should, strictly speaking, be the same 
for two codes meeting at a given interface plane. 

• X2: Axial location of the interface plane downstream of the acoustic element relative to the origin 
(non -dimensionalized by DFAN) as specified by the user. 

• DIA2: Outer duct diameter of the interface plane downstream of the acoustic element 
(non— dimensionalized by DFAN). 

• MX2: Mass averaged axial Mach number at the interface plane downstream of the acoustic 
element. 

• MTIP2: tip rotational Mach number at the interface plane downstream of the acoustic element. 

• PS2: Mass averaged static pressure at the interface plane downstream of the acoustic element 
(non— dimensionalized by PINF). 

• TS2: Mass averaged static temperature at the interface plane downstream of the acoustic element 
(non— dimensionalized by TINF). 

Card 12: NDAT: Number of Streamlines where geometry is specified 

Card 13: DLA(I), 1=1, NDAT: Streamline diameter (non— dimensionalized by DFAN) 

Card 14: SSOLIDITY(I), 1=1, NDAT: Stator Solidity (stator chord/stator pitch) 

Card 15: SALPHCH(I), 1=1, NDAT: Stator Stagger Angle relative to the circumferential direction 
(degrees). This is normally considered to be a number between 0 and 110 degrees. 
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Card 16: SSOC(I), I=1,NDAT: Upstream rotor trailing edge to stator leading edge aerodynamic 
Spacing to upstream rotor aerodynamic chord for each streamline (Set equal to zero if there 
is no upstream rotor). 

1.3.4.2 Scattering Coefficients (see Section 1.3.6) 

Card 17: Write: FORMAT( IX, ’Scattering Matrix and Source Vector’) 

Card 18: Write: FORMAT( IX, ’Scattering Matrix’) 


Card 19: Scattering coefficients as follows 


C 

C 

c 

c 

c 

c 

c 


WRITE SCATTERING MATRIX INFORMATION S 12 , S13, S14, AT T HE 
INTERFACE PLANE UPSTREAM OF THIS REGION. THESE SCATTERING 
COEFFICIENTS CORRESPOND TO DOWNSTREAM PROPAGATING WAVES 
ENTERING THE UPSTREAM INTERFACE PLANE AND BEING REFLECTED 
BY THE ELEMENT BACK OUT AS UPSTREAM WAVES. 

DO 70 IW = NUPSTRE + 1 ,NWAVETYPE 
DO 80 IWP = 1, NUPSTRE 


& 

MU = 0, MUMAX(N,K,IW)), 


& 

MUP = 0, MUMAX(N,KP,IWP )), 


& 

K = KMIN(N,IW), KMAX(N,IW)), 

k 

& 

KP = KMIN (N,I WP),KMAX(N ,IWP )) , 

k' 

& 

N =1, NHT) 

<r~ n 

80 

CONTINUE 


70 

CONTINUE 



C 

C 

C 

C 

C 

C 

C 


WRITE SCATTERING MATRIX INFORMATION, S n ,AT THE INTERFACE 
PLANE UPSTREAM OF THIS REGION. THESE SCATTERING COEFFICIENTS 
CORRESPOND TO UPSTREAM PROPAGATING WAVES AT THE 
DOWNSTREAM INTERFACE BEING TRANSMITTED THROUGH THE 
ACOUSTIC ELEMENT TO THE UPSTREAM INTERFACE 

DO 90 IW = 1, NUPSTRE 
DO 100 IWP = 1, NUPSTRE 
WRITE(IUNIT,999)((((( Sj^N, K, KP, MU, MUP), 

& MU = 0, MUMAX(N ,K,IW)), 

& MUP = 0, MUMAX(N,KP,IWP)), 

& K = KMIN(N,IW), KMAX(N,IW)), 

& KP = KMIN(N,IWP),KMAX(N,IWP)), 

& N =1, NHT) 

100 CONTINUE 
90 CONTINUE 


<-* 

n 
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c 

C WRITE S22, S32, S42, S23, S33, S 43 , S 24 , S34, S44. THESE COEFFICIENTS 
C CORRESPOND TO DOWNSTREAM PROPAGATING WAVES INPUT INTO THE 

C REGION AT THE UPSTREAM INTERFACE PLANE WHICH ARE 

C TRANSMITTED THROUGH THE ELEMENT TO THE DOWNSTREAM 

C INTERFACE PLANE 
C 


DO 90 IW = NUPSTRE+ 1,NWAVETYPE 
DO 100 IWP = NUPSTRE + 1 ,NWAVETYPE 

WRITE(IUNIT,999)((((( SgJ^N, K, KP, MU, MUP), 
& MU =0, MUMAX(N,K,IW)), 

& MUP = 0, MUMAX(N,KP,IWP)), 

& K = KMIN(N,IW), KMAX(N,IW)), 

& KP = KMIN(N,IWP),KMAX(N,IWP)), 

& N = 1, NHT) 

100 CONTINUE 
90 CONTINUE 


+-k 

<-n 


C 

C WRITE SCATTERING MATRIX INFORMATION, S 2 i , S 31 , S 4} , AT THE 
C INTERFACE PLANE DOWNSTREAM OF THIS REGION. THESE SCATTERING 
C COEFFICIENTS CORRESPOND TO UPSTREAM PROPAGATING WAVES 
C REFLECTED BY THE ACOUSTIC ELEMENT AS DOWNSTREAM 
C PROPAGATING WAVES. 

C 


DO 110 IW = 1, NUPSTRE 
DO 120 IWP = NUPSTRE + 1 ,NWAVETYPE 


WRITE(IUNIT,999)(( ((( S#P„(N, K, KP, MU, MUP), 
& MU =0, MUMAX(N,K,IW)), 

& MUP = 0, MUMAX(N,KP,IWP)), 

& K = KMIN(N,IW), KMAX(N,IW)), 

& KP = KMIN(N,IWP),KMAX(N,rWP)), 

& N = 1, NHT) 

120 CONTINUE 
110 CONTINUE 


k 

< h - n 
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I.3.4.3 Source Vector (see Section 1.3.6) 


Card 20: IFORCE = 0 No Source Vector is Input 
= 1 Source Vector is input 

IF IFORCE = 1: 

Card 21: Write: FORMAT(lX,’Source Vector’) 

Card 22: Upstream propagating pressure mode amplitudes normalized by PINF 

DO 10 IW = 1,NUPSTRE 
READ(IUNIT,999) ((( B%( MU,K,N), 

& MU = 0, MUMAX(N ,K,IW)), 

& K = KMIN(N,IW),KMAX(N,IW)), 

& N =1, NHT) <"« 

10 CONTINUE 

Card 23: Downstream propagating pressure mode amplitudes normalized by PINF 

DO 10 IW = NUPSTRE + 1,NPRESDIR 
READ(IUNIT,999) ((( J?J(MU,K,N), 

& MU = 0, MUMAX(N,K,IW)), <-f* 

& K = KMIN(N,IW),KMAX(N,IW)), <-k 

& N =1, NHT) «-» 

10 CONTINUE 

Card 24: Downstream propagating vorticity waves normalized by c^. 

DO 10 IW = NPRESDIR + 1, NWAVETYPE 
READ(IUNIT,999) ((( R^(MU,K,N), 

& MU = 0, MUMAX(N,K,IW)), 

& K = KMIN(N,IW),KMAX(N,IW)), 

& N =1, NHT) «"» 

10 CONTINUE 

END IF 
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1.3.5 Aft Radiation Code: 20 < ICODE < 30 
I.3.5.1 Header Cards 


Information in Cards 11 through 16 are only used to confirm for the user which Acoustic Properties 
File the user has chosen. These parameters are not used in code calculations. While parameters in 
Card 10 are not presently required for code calculations, this information expected to be required for 
code calculations in the future. 

Card 10: 

• SIGMAR: Hub to tip ratio of the duct at XINTER. This number should, strictly speaking, be the 
same for two codes meeting at a given interface plane. 

• XINTER: Axial location of the interface plane relative to the origin (non-dimensionalized by 
DFAN) as specified by the user. 

• MINF: Far-field axial Mach number (far from the nacelle, positive downstream) 

• MDUCT: Mass averaged axial Mach number at the interface plane (positive downstream) 

Card 11: NLINO: Number of acoustically Lined Elements (outer wall). Set this value equal to zero 
if there is no outer wall liner. 

If ( NLINO .GT. 0 ) THEN 

Card 12: 

• MBEGO: Outer wall element number for start of liner 

• XBEGO: Axial location where outer liner begins relative to the origin (non— dimensionalized by 
DFAN) as specified by the user. 

• XENDO: Axial Location where outer liner ends relative to the origin (non-dimensionalized by 
DFAN) as specified by the user. 

Card 13: WRITE(IUNIT,999) ZIMPO (I), I = MBEGO, NLINO + MBEGO: Liner impedance 

(Resistance, Reactance) for each element on outer wall non-dimensionalized by local qc (i.e. 
density x acoustic speed). 

ENDIF 

Card 14: NLINI: Number of acoustically lined elements (inner wall). Set this value equal to zero if 
there is no inner wall liner. 

If ( NLINI .GT. 0 ) THEN 

Card 15: 

• MBEGI: Inner wall element number for start of liner 

• XBEGI: Axial location where inner liner begins relative to the origin (non-dimensionalized by 
DFAN) as specified by the user. 

• XENDI: Axial Location where inner liner ends relative to the origin (non-dimensionalized by 
DFAN) as specified by the user. 


42 



Card 16: WRITE(IUNIT,999) ZIMPI(I),I=MBEGI,NLINI+MBEGI: Liner impedance 

(Resistance, Reactance) for each element on the inner wall non— dimensionalized by local gc 
(i.e. density x acoustic speed). 

END1F 


I.3.5.2 Scattering Coefficients (see Section 1.3.6) 

Card 17: Write: FORMAT(lX, ’Scattering Matrix and Source Vector’) 
Card 18: Write: F0RMAT(1X, ’Scattering Matrix’) 


Card 19: Scattering coefficients as follows 
C 

C WRITE SCATTERING MATRIX INFORMATION, S }2 , AT THE INTERFACE 
C PLANE UPSTREAM OF THIS REGION. THESE SCATTERING 
C COEFFICIENTS CORRESPOND TO DOWNSTREAM PROPAGATING WAVES 
C SCATTERED INTO UPSTREAM PROPAGATING WAVES. NOTE THAT 
C Sji& S 14 ARE NOT WRITTEN SINCE THEY ARE VORTICITY WAVE BASED. 
C VORTICITY WAVES ARE NOT MODELED AS INPUT TO THE AFT 
C RADIATION CODES. THUS S 13 AND S 14 ARE ZERO. 

C 

DO 10 IW = NUPSTRE+ 1,NPRESDIR 


c 

DO 20 IWP = 1, NUPSTRE 
WRITE(IUNIT,999)(((( Sj^N,K,MU,MUP), 


& 

MUP = 0, MUMAX(N,K,IWP)), 

& 

MU = 0, MUMAX(N,K,IW)), 


& 

K = KMIN(N,IW),KMAX(N,IW)), 

<r- k 

& 

p 

N = 1, NHT) 

n 

L 

20 

CONTINUE 


10 

CONTINUE 


I.3.5.3 

Source Vector (see Section 1.3.6) 



Card 20: IFORCE = 0 No Source Vector is Input 
= 1 Source Vector is input 


IF IFORCE = 1: 

Card 21: Write: FORMAT(lX, ’Source Vector’) 

Card 22: Upstream propagating pressure mode amplitudes normalized by PINF 

DO 10 IW = 1,NUPSTRE 
READ(IUNIT,999) ((( B^MU,K,N), 

& MU =0, MUMAX(N,K,IW)), <-/< 

& K = KMIN(N,IW),KMAX(N,IW)), <-k 

& N =1, NHT) <-* 

10 CONTINUE 
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Card 23: Downstream propagating pressure mode amplitudes normalized by PINF 

DO 10 IW = NUPSTRE + 1, NPRESDIR 
READ (IUNIT,999) ((( R^(MU,K,N), 

& MU =0, MUMAX(N,K,IW)), 

& K = KMIN(N,IW),KMAX(N,IW)), 

& N = 1, NHT) <r-n 

10 CONTINUE 

Card 24: Downstream propagating vorticity waves normalized by c^. 

DO 10 IW = NPRESDIR + 1, NWAVETYPE 
READ(IUNIT,999) ((( 5^(MU,K,N), 

& MU =0, MUMAX(N,K,IW)), 

& K = KMIN(N,rW),KMAX(N,IW)), <-k 

& N = 1, NHT) <-n 

10 CONTINUE 

ENDEF 

I.3.5.4 Far-field Directivities 

Card 25: Write: FORMAT(lX,’Farfield Directivities’) 

• NANGLE: Number of polar angles ( 6 in Figure 5) in the directivity pattern. The code will 
interpolate these angles to form an array every 1° from 0°< ANGLE < 180°. In the present code, 
any angle which is within 25° of the forward most computational boundary will be discarded. This 
is consistent with the baffle boundary condition which presently is found in the Eversman Aft 
Radiation code. 

• DARRAY: Radial distance (R in Figure 5) from the origin non-dimensionalized by DEAN 
(origin specified by the user) of the polar directivity array. 

Card 26: Far— field angle in degrees (0° is in front of the inlet) 

• WRITE(IUNIT, * ) ( ANGLE(IANGLE), IANGLE = 1, NANGLE) 

Card 27: Write: FORMAT( IX, ’COMPLEX PRESSURE DIRECnVITIES’) 

Card 28: Complex pressure directivities (real and imaginary parts of the directivities for each mode): 

DO 10 IW = NUPSTRE + 1, NPRESDIR 
C 

WRITE(IUNIT,999) (((( DIRECTION, K,MU, IANGLE), 

IANGLE = 1, NANGLE), 

MU = 0, MUMAX(N,K,IW)), 

K = KMIN(N,IW), KMAX(NJW)) 

N =1, NHT) 

10 CONTINUE 
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1.3.6 Format for Scattering Coefficient and Source Vector Output 

Section 3.2 discusses the conventions for the scattering coefficients and source vector information in 
CUP3D. Further discussion of the “Kinematics of Mode Scattering” may also be found in Reference 
1. The discussion in this section will center on how scattering coefficients and source vectors are 
non-dimensionalized and represented in CUP3D. 

A two sided Fourier transform is utilized where the complex scattering coefficient and source vector 
amplitudes represent values of the coefficients in this series. Equations (20) through (25) in Section 
3.2 give further details. Only the positive harmonics of the source vectors and scattering coefficients 
are input. Because the code works with tone noise, the negative harmonics are assumed by the code 
to be the complex conjugates of the positive harmonics with the exception of rotor scattering. In the 
rotor case, positive harmonics scatter into negative harmonics thus requiring scattering coefficients 

to be input. 

All scattering coefficients, source vectors and far— field directivities are read in with format given by. 
999 FORMAT(10E13.6) 

where this format will include five coefficients (i.e. real and imaginary parts) per line. 

The scattering coefficients, S, are effectively the reflection and transmission coefficients for each 
acoustic element. The scattering coefficients take on the form S P J W where the scattering coefficient 
subscripts denote the wave type, W, being scattered into the wave type, W’. The superscripts denote 
interface plane, P, scattering into the interface plane, P . 

All source vector information should be non— dimensionalized by: 

• The ambient pressure, PINF for pressures waves, 

. The ambient acoustic speed, Coo = frJJZ where T M = TINF, for vorticity waves. 

TINF and PINF are defined as the static temperature and pressure far from the nacelle. 

In the present code, the following wave types are in use: 

W = 1 Upstream going pressure waves 
W=2 Downstream going pressure waves 
W = 3 “no-w” vorticity wave (downstream only) 

W = 4 “no-u” vorticity wave (downstream only) 

In Sections 1.3.2 to 1.3.5, FORTRAN comment statements are shown which refer only to these four 
wave types. However, the number of wave types can be changed by changing NUPSTRE, 
NPRESDIR, and/or NWAVETYPE in the parameter statements in the CUP3D code. 

For scattering coefficient and source vector input formats discussed in Sections 1.3.2 to 1.3.5, the 
superscripts (P and P’) are given by: 

U = interface plane upstream of acoustic element = IUPSTR(IREGION) in CUP3D System File 
D = interface plane downstream of acoustic element = IDWNST(IREGION) in CUP3D System File 
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For scattering coefficient and source vector input formats discussed in Sections 1.3.2 to 1.3.5 the 
subscripts (W and W’) are given by: 

W = IW = Wave type of the wave travelling into the acoustic element 
i.e. right hand side of Equation (15) 

W’ = IWP = Wave type for the wave leaving the acoustic element (i.e. the scattered wave) 
i.e. left hand side of Equation (15) 

NUPSTRE = Number of upstream travelling pressure wave types (= 1 in the present code) 
NPRESDIR = Total number of pressure wave types (= 2 in the present code) 

NWAVETYPE = Total number of wave types: pressure + vorticity (= 4 in the present code) 
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1.3. 7 Examples of Acoustic Properties Files 

Example 1: Inlet Radiation Code (same style as aft radiation code) 
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APPENDIX II: POWER CALCULATION IN CUP3D COUPLED 

CASCADE CODE 


The objective of this appendix is to develop the equations being used to calculate power in CUP3D. 
The sound power level is calculated at interface planes just upstream and just downstream of the noise 
source. The following assumptions are made at those interface planes: 

• Locally constant area annular duct, 

• Locally uniform axial flow 

As a result, for the sound power calculation, the standard Tyler/Sofrin modes (Reference 13) are used. 
A more general formulation may be developed in the future. 

From Equation (4), the perturbation pressure is assumed to be of the form: 


oo oo 


n = — oo k = — oo ^ = 0 

Perturbation velocity may be related to perturbation pressure by: 


(III) 


(II.2) 


n- ~ oo k = — oo [x =0 


R Wnkp ma y t ^ en be found using the small perturbation, unsteady axial momentum equation for a 
uniform axial mean flow (consistent with Equation C.ll in Reference 11): 


% + U = 0 

dt ? fa PpdX 

Substituting Equations (II.l) and (II.2) into Equation (11.3) we get: 

!C * ( ~ + UpyWnkfi)- 

Y 

Multiply through by - and solve for R^ nlcM to get: 


\ dp 

'^WnkfL 


- iP?- V P 

1 pp'Wnkfi 


(II.3) 


(II.4) 


Where: 



P oo 

C oo P pC p 


R Wnkp 



(H5) 


(II.6) 


Note that for propagating waves, R^ nkfi is real. This expression is consistent with Equation C.14 in 
Reference 11 except that Reference 11 non-dimensionalizes differently. 
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The equation for time dependent intensity in a uniform axial flow is used i.e. Equation C.4 in 
Reference 11 or Equation 101 in Reference 1: 


if{x,r,<p,t ) = (i + Miyy w 


p^ipw) + Pj£pMx P {uw) 


(II.7) 


Upstream acoustic power will be calculated at the interface plane, P = 1 using upstream propagating 
acoustic waves (for a uniform axial flow there is one upstream going acoustic pressure wave type: W 
= i) Downstream acoustic power needs to be calculated at the last interface plane, P — N p , using 
downstream acoustic pressure waves (for a uniform axial flow there is one downstream going acoustic 
pressure wave type: W = 2). Acoustic power is expressed by. 


n P = 


r d 1 2jc 

Ijfdcpdr = r 2 d j j Tprd<pdr 

r h 0 ° 0 


(II 8) 


where f is the dimensional radius and I P is the time average intensity of I P . To get I P , realize that only 
cuton modes contribute to the sound power level. In Equations (H.l) and (II.2) the perturbation 

pressure, pL and perturbation velocity, u p w are real. This means that p p w = p p w and u p w = . 

If Nh is given as the number of BPF harmonics being calculated and only cuton modes are included 
in the summation, then substitution of Equations (II.l) and (II.2) into Equation (II.7) leads to: 


I P (X,r,<p,t) = 




(II.9) 


---N h k p n' = -N h k' p‘ 


where: 


\nkpm'k'fi' 


(l + M^jpooCooR^,^, + + PFCpMxpC^WnlcfPwnlc', 


(11.10) 


To get the time average intensity at interface plane, P, note that: 


, —i(n -ri)BQt _ ^ 


(ini) 


so that Equation (II.9) becomes: 


N, 


Ip[Xp>r,(f>) = XXX lL A Wnk P A Wnk' P 'hnk P k' p y r w'P 

n--N h k pi k f n' 


m pi 


te i(Y-y'){x-x P ) e i(k-k’)y<t>^ lU2 ) 


where the subscripts on I x in Equation (II. 10) has been contracted in the obvious way. Substitute this 
equation into Equation (II. 8): 


1 2jt 

1111 f f (il- 13 ) 

KJ t, u' J J 
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The integral over <f> yields a Kronecker delta: 


2 71 

| e -i(,k-k')v<t> d( p _ 2 ^. (11.14) 

o 


so that: 

N h r 

= t (ins) 

n=-N h k fi n‘ J 

The orthogonality condition for the radial eigenmodes is: 

l 

j 'Pn^fp^rdr = d^,C mfl (11.16) 

O 

Thus when fi ' = fi: 

l 

^ TpmfjTdr — C m<u (11.17) 

a 

This integral is evaluated in Appendix III to determine the value of C mil . In addition, since negative 
harmonics are equal to positive harmonics, use positive harmonics only and multiply Equation (11.15) 
by 2. The Kronecker delta in Equation (11.16) enables the p' summation and eliminates the final 
exponential in Equation (11.15) since all of the implied subscripts on y and y’ are now equal. Therefore 
Equation (11.15) becomes: 


N k 


n P 4jIr d X X ^ A W„kv A W„kJlnkti C ' 

n = \ k fi 

Equation (II. 10) can now be substituted back into Equation (11.18): 

HP = X X Y^WnkAwn 


mfi 


n = \ k fi 


WnkfTWnkfi 

M 


(l + Mx p )p <x>c coR^ nkfl + pp/pP 2 * + P F £ pMx p clc,Rw nkfl RwnkM 


" mfi 


(11.18) 


(11.19) 


Non — dimensionalize power by p « c » (2r d ) so that: 




”p = *22 2<^ A 


n = l k fi 


Wnkfi Wnkfi 


(* + KK2 + Cm (H-20) 
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Substitute Equation (II.5) into Equation (11.20) to get: 


Nu 


tip = Z Z Z A Wnkp A Wnkp 


n - 1 k fi 

(l + Ml^Rwnlqi + M *r + M x^V n k t fiw n kp jCw 


( 11 . 21 ) 


Since the y’s are 
(11.21) becomes: 


real for propagating modes, Equation (II.6) shows that Rw = Rw and Equation 


N, 


h |2 

l P ^ PpCfC 00 /-J 2 ^\ Wnkfx] 

n = l k ii 


Tip = 71- 


(l + M 2 Xp )R P Wnkfi + M Xp + Mx^Rw^J 


'mpt 


(11.22) 


Since c z = — = y a R a T where P = mean static pressure, Equation (11.22) can be written: 






n — \ k fi 


P 

Wnkp 


(l + Ml^jR-Wnkfi + M *P + M x P (kw„k^ 


C m u (11.23) 


This equation gives the sound power for blade passing frequency harmonics combined up to harmonic, 
N h . The limits for k and/* are set by the number of cuton modes. 


This equation can also be expressed as the sum of modal power levels, TI Pnkfl , where: 


N, 


n F ~ 

n = l k fi 


(11.24) 


and: 


tv r 00 

n Pnkti = ^YaPp 



r “K^I 


(l + M z p )Rw nkfl + M Xp + M Xp ^R P Wnk ^ 


'mu 


(11.25) 


This equation is equivalent to Equations (33) and (34). Rw^ is found in Equation (II. 6) consistent 
with Equation (35), and Cmp is found in Appendix III. 

Thus, sound power level for a particular harmonic is: 

= < IU6) 

k p 

where the summation is over propagating modes. 
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APPENDIX III: NORMALIZATION OF THE UNIFORM FLOW 

DUCT MODES 

The purpose of this appendix is to evaluate the integral in Equation (11.17) analytically. Start with the 
expression for Bessel’s equation for the radial direction. This equation is developed by applying the 
method of separation of variables to the converted wave equation in a mean axial flow. The equation 
related to the radial direction corresponds to Equation 4.4 in Reference 11. 

+ = 0 (nu) 

where k is the separation constant and r is the dimensional radius. It will later be shown to also be 
the mode eigenvalue. 


With the notation ip' = 


_ dip. 


dr' 


np" + ip’ + | k z r - %^ip = 0 

The boundary conditions for this equation for a hardwalled annular duct are: 


(IH*2) 


Multiply Equation (III.2) by 2np'(icr): 


1p'{KT d ) = 0 
ip'[Kr h ) = 0 


(III.3) 

(III-4) 


2 ?ip"ip' + 2Hip') 2 + 2(?c 2 r 2 - m 2 'jipip' = 0 (III.5) 

Noting that ) 2 J = 7F 1 ip''ip' + 2f(ip') 2 and (ip 2 '} = 2ipip' Equation (III.5) becomes: 

(/^') 2 ] = - (jcV - m 2 )(v> 2 )' (III.6) 

Integrate both sides of this equation across the duct: 

r d r d 

j |^('V'') 2 j dr = - j (k 2 ? 2 - m 2 )(v> 2 ) dr (HL7) 

r h r h 

or: 

r d 

j = - j (k 2 * 2 - m 2 )(ip 2 ) dr (III-8) 

r h 

By applying boundary conditions Equations (III.3) and (III.4) to Equation (III.8), the left hand side 
of this equation becomes zero and the ip on the right hand side becomes the eigenfunction ipm^mf/), 
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within a constant. However, we do not apply the subscripts until the normalization is defined below. 
Therefore: 


0 = - 




(HI*9) 


Integrating the right hand side by parts and rearranging: 


1 

2 





or: 


(III. 10) 


1 

2 





(HI 11) 


Factor out r^, normalize tc by and normalize r by r^. Also normalize ip by VWu where V-'max is the 
maximum value of each mode eigenfunction (Vw = 1 in CUP3D). Thus: = Kr d , 


a = 



hub/tip ratio of the duct, and r 


— is the non-dimensional radius. Therefore: 

r d 


1 

2 





(III.12) 


The integral on the right hand side of this equation is the same as the left hand side of Equation (11.17). 
Thus C mft is: 




1 



(III.13) 


This equation applies to a locally constant area hardwall annular duct with locally uniform axial flow. 

are the mode eigenvalues and are the mode eigenfunctions. Note that to calculate power 
only the mode eigenfunctions at the hub and tip locations are required. 
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APPENDIX IV: CALCULATING THE TOTAL ACOUSTIC 

FAR-FIELD PRESSURE 


The purpose of this section is to derive the equations for the far-field directivity using the state vector 
solution (i.e. mode amplitudes) from CUP3D. CUP3D assumes that the system is formulated in the 
frequency domain where acoustic pressure is represented by the superposition of ( m , p) modes for 
each blade passing frequency harmonic. In polar coordinates of Figure 5 acoustic pressure can be 
expressed by: 


p(R,e,<p,t) = x X Y,Pnk fi ( R ^y (m,p ~ nBat) (iv. i) 

n = -oo k- - oo ^ = 0 

where p ^ are the complex pressure mode amplitude coefficients obtained at the desired far-field 
points from a radiation code. Since negative harmonics are the complex conjugate of positive 
harmonics (i.e. =p^): 

{ 00 00 OC | 

X X XPn^ R ’ d y (m4> ~ nBQt) \ (iv.2) 

tt = lfc=-oo^=0 J 

Equation (IV. 1) as applied to the far-field can also be expressed as the combination of pressure waves 
propagating from the inlet and nozzle. Thus: 

= PlnkM ^’ + PAnk^R’ *) (IV.3) 

where p Inkfl (R,Q) is the inlet complex pressure directivity for each (m, p) mode, and p Ankf jR,6) is 
the aft complex pressure directivity for each (m,p) mode. These complex pressure directivities can 
also be expressed as a superposition of state vector amplitudes, A multiplied by directivity shape 

functions, D , ^ , at interface planes just upstream and just downstream of the noise source, for each 
radial mode order. 

For inlet noise the interface plane, P, where the inlet is connected to the noise source is used (e.g. in 
Figure 7 ,P = 1) to obtain the upstream propagating pressure waves from the state vector, A^ nkM and 
the mode inlet directivity shape functions, D IfVnkft (R, 0). The CUP3D code is configured to permit 
multiple upstream pressure wave types so that these pressures are added together to obtain the 
complex pressure directivity for each (m, p) mode, p Inkfi (R,6). The number of forward propagating 
pressure wave types is given as N UPS jr ■ Thus: 

^UPSTR 

PlnkfA R ’Q} ~ X r nk/t( R ’@) (FV.4) 

^=1 

Note that for all modal representations developed to date (e.g. Reference 14) only one type of forward 
propagating mode has been found (i.e. N UPSTR = 1). However, CUP3D uses Equation (IV.4) as shown 
to permit more upstream pressure wave types to be added in the future, if necessary. 

For aft noise the interface plane where the aft radiation code is connected to the noise source is used 
(e.g. in Figure 7, P = 3) to obtain the downstream propagating pressure waves from the state vector, 
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Aw nkfi and the mode aft directivity shape functions, D AWnkfi (R,0). The CUP3D code is configured 
to permit multiple downstream pressure wave types so that these pressures are added together to 
obtain the complex pressure directivity for each (m,/u) mode, Pxnkfi{ R ’^)- The num h er °f propagating 
pressure wave types is assumed below to be Npresdir so that the aft pressure directivity for each (m,fx) 
mode becomes: 

N PRESDIR 

PAnkfS- R = X A WnkfPAWnkM^ R ’ e ^ ^ IV5 ' > 

W= ^UPSTR + 1 

For the modes presently in TFaNS and discussed in Section 3.2, Npresdir = 2 - However, CUP3D uses 
Equation (IV.5) as shown to permit more downstream pressure wave types to be added in the future, 

if necessary. 

Thus, to calculate the far-field directivity the following items are needed for each pressure wave type: 

1. The solution to the state vector, A p Wnkfi , at the interface plane where a radiation acoustic element 
is connected to the noise source, inlet and aft. 

2. Directivity shapes for each propagating pressure wave type from an inlet radiation calculation, 
D JWnkfl (R,0) and an aft radiation calculation, D AWnkfl ( R ,9)- These are determined by inputting 

unit mode amplitudes into a radiation code and then calculating the far— field directivities. 

Also, CUP3D calculates the following far-field pressure directivities: 

• Total far— field sound pressure level (SPL) directivity (inlet and aft combined) 
for each harmonic 

• Inlet SPL directivity for each harmonic 

• Aft SPL directivity for each harmonic 

• Inlet, aft and total SPL directivities for each circumferential mode order at each 
harmonic (summed over all radial mode orders) 

• Inlet, aft and total SPL directivities for each circumferential and radial mode 
order at each harmonic 

For this derivation we will concentrate on total directivities which can be represented as the 
combination of the inlet and aft directivities. Then more simplified forms of the equations will be 
derived. Given Equations (IV.l), (IV.3), (IV.4), and (IV.5), the acoustic pressure in the far-field at 
a given directivity angle in CUP3D is: 

P{R,e,4>,t)= X X X Pink fi ^ e y im(p ~ nBQt) + Pm^ r > 0)e^- nBQt) (iv.6) 

n = — co k= - 

Or given that m = nB - kV: 

KMM- 1 1 <- wr > 

n — — <x> k= — oo ^ =0 

Let O nk (<t>,t ) = [nB - kV)(p - nBQt so that Equation (IV.7) becomes: 

P(R,6,4>,t) = X X Xh^ + PAnkpY 6 ^ < IV8) 

— m L- — — so r; 
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This equation now needs to be processed to obtain the root mean square acoustic pressure as follows: 


P(R,6) = Jp 2 (R,6,<P,t) (IV.9) 

The notation imphes the following operations: first square P, then time and circumferentially average 
P 2 , and then take the square root. A circumferential average is being calculated because we rarely 
know the engine azimuthal angle relative to the location of engine circumferential origin. Thus: 

P 2 = PP'= X X X X X 2>*-, + (IV.10) 

n = — oc fcs* — oo Jc'« — » ^'*0 

Substituting back in for 0 n k(<P>O and rearranging the exponent gives: 

p2 = X ii X X )^i (IV11) 

n* - oc - go ^ «0 n’ = — oo k' * oo ^'*0 

The magnitudes in the above equation can also be expressed as the summation over all /z’s. Thus 
rearrange this equation to get: 


' 2 = X X X X < 

rtse-oofc=-ao/j’*:-ooA:'--» 


00 

^ P Anku 




yy + PAn'k'p' 


„i\\(n-n')B-(k—k r )V)fi-(n-n' )BQt 


i(rv.i2) 


Now average over time. It can be seen that when n^n\ then the acoustic pressure is sinusoidal and 

acoustic pressure over time averages to zero. Thus we get e i ^ n ~ n ' iBQt = d^, and Equation (IV.12) 
becomes: 


p\R,e,M- £ £ £ 


n=-«fca= — 00 it'* - « 


M =0 


w 00 

^ Inku Adnty] ' /* ' [P/aA'^ f 

J 


e -i|(* -*')»> (IV.13) 


Let: 


Pn* - i 


00 

X. [-P /nfat ’ 

U=° J 


(IV.14) 


and: 


PnJt’ ~ “j X [P Ink 'P' + PAnky\ 

U'=0 


So that Equation (IV.13) becomes: 


(IV. 15) 


P 2 (R,d,<p,t ) = X X 2 PnkP'nk'e- i{(k - k ' )V} + (IV.16) 

n~ — <x> k— — ocjt'= — oo 

This can also be expressed as a one sided series on n where: 

PHR, 0,0,0 = X X X (IV. 17) 

n = lk = — »<:’ = — oo 
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where we have dropped the n = 0 term. In the second term note that: p and 
P - n k = Pn-k SO that. 


00 00 00 


P\R,e,M = 2^ X X ArfP** ' 


i{(k-k')V)<p 


(IV. 18) 


n = lk = - » A' = - ® 


Equation (IV.18) shows that the root mean square acoustic field is spatially periodic in the <p 
(circumferential) direction with the same number of lobes as vanes. Substituting back Equations 
(IV. 14) and (IV. 15): 

7 ) -2jr X X fe- c™ 19 ) 

„«i ^«o J =° 

For our purposes we want the mean square pressure directivity for each harmonic separately so that: 


P 2 n (R,0,4>,t ) = 2 


x x (iv20) 


The square root: 


P 2 n (R,e,<t>,t ) = / 2 


t*-oo *'«-« J |/~° J 


is the general equation for the root mean square pressure directivity at an arbitrary circumferential 
angle. This equation is useful if the engine azimuthal angle is known relative to a known 
circumferential origin. In general the circumferential origin and the azimuthal angle may be known, 
but the effect of phase shifting in the noise measurement equipment is not known accurately. Thus, 
it is desirable to average Equation (IV.20) circumferentially. In Equation (IV.2 0), when k # k’, then 

the pressure is sinusoidal and averages to zero circumferentially. Therefore e~ l(k ~ k )V<t> = <5**' so 
that: 


P 2 n {R,6,<P,t ) =P n P' n = 


C f 

Wn = 2 X < ^[Plnkp+PAnkfi] M X W + P**A ' ( IV22) 

k=-<x> ji = 0 * J** m Q - 


Since the summation over p is the complex conjugate of the summation over p’: 


Pl{R, 6,<p,t) = X X & [ p MM + p A*\ 

k= - » [u =0 


(IV.23) 


To get Equation (IV.9) take the square root of Equation (IV.23) and rearrange to get Equation (41) 
in Section 3.4.2: 


P n (R,6)= / X 

\j k= — “^ = 0 fi= 0 
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This gives us the root mean square of the total harmonic pressure pattern averaged circumferentially. 

Note that pressure as input to CUP3D is normalized by a reference pressure, p^. Since 
to get far— field SPL Levels: 


SPLntotal 


= 20 log ^ — — 
\ F ref 


) 


where P re f = 20pPa. 

Simplified versions of Equation (IV.24) are discussed in Section 3.4.2. 


(IV.25) 


64 



APPENDIX V: MODIFYING EVERSMAN INLET AND AFT 

RADIATION CODES TO INCORPORATE THEM 
INTO TFaNS 


The Eversman Inlet and Aft Radiation Codes have been incorporated into TFaNS to calculate far 
field noise and provide internal reflections from inlets and nozzles. The purpose of this appendix is 
to discuss the modifications to these codes and to derive the equations used in the modified codes. 
Modifications to the radiation codes include: 


• Capability to run multiple modes and harmonics with unit input. 

• Dimensional input format based on standard geometry and performance param- 
eters. 

• An Acoustic Properties File output for use by CUP3D. 

In the sections which follow, these modifications will be further discussed. 


V.1 MULTIPLE MODE RUNNING & CHANGES TO THE INPUT FORMAT 


The capability to run multiple modes and harmonics was added by creating a mode search routine. 
This routine decides which mode and frequencies to run through a radiation module and then runs 
the module with unit input for each mode. 


As part of the creation of TFaNS versions of the inlet and aft radiation modules, the input file format 
was changed to be more consistent with standard design parameters. Two parameters are now calcu- 
lated instead of being input: The reduced frequency, ETAR (*) and the circumferential mode order, 
m. The reduced frequency is defined by: 


V r = 


C 00 


(V.l) 


Since a> = nBQ c » 


= JvaRaTo,, and Q = 


2itN u 

~60~ : 


n r - nB 


360 JyJR^T „ 


(V.2) 


where: 

Ni = fan rotational speed (rpm) 
r d = fan duct radius ( inches ) 

R a -Ideal gas constant for air = 1 716.5 ft 2 /sec R 

y a = ratio of specific heats = 1.4 for air 

In an adiabatic flow the total temperature and static temperature can be related to the axial Mach 
number far from the nacelle as: 


l£i= i + *Ll 
T x 2 


-ML 


Therefore substituting Equation (V.3) into Equation (V.2): 


Tj r = nB 


/ \ 1/2 

i + *yV.) 

360 ly a R a T 0i \ 2 / 


(V.3) 


(V.4) 
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A common design parameter is the corrected fan rotational speed, N u : 


Substituting Equation (V.5) into Equation (V.4): 

Vr = nBM tipc 

where: 


flhL 

V 518.67 


(V.5) 


(V6) 


M 


tipc 


?Mi J 1 + vjlz± mL 


1/2 


360v/y a i? fl (518.67) V 


(V.7) 


Equations (V.6) and (V.7) are used by both the inlet and aft radiation codes. N lc , r d and B are input 
by the user. M x „ is the Mach number far from the nacelle and is available from the radiation code 
mesh generation output file. 

Also calculated by the TFaNS radiation modules is the circumferential mode order, m = nB —kV. This 
is accomplished by having the user input the number of stator vanes, V into the code. The mode search 
routine then determines the modes it needs to calculate by varying k and fx for a user specified number 
of BPF harmonics, n. 


V.2 MODIFYING RADIATION CODES TO OUTPUT ACOUSTIC 
PROPERTIES FILES 


To correctly output information to an Acoustic Properties File, two major items must be output from 
the radiation codes: 

• Scattering Coefficients 

• Far— field directivities 

A method for creating this output is discussed below. 

V.2.1 Scattering Coefficients 

Both the inlet and aft radiation codes are formulated from the same theory using the same basis for 
a coordinate system and exponential notation (see Reference 3). The inlet and aft radiation codes 
are formulated in terms of the perturbation velocity potential. Acoustic pressure can then be calcu- 
lated from this velocity potential. For a single acoustic mode at the radiation code input plane, the 
perturbation velocity potential, <p, is given by Equation 34, Reference 3: 

/ Wnkfi (X,r,<p,t) ~ (V.8) 

where the notation has been made consistent with this report except for kx which is the axial wavenum- 
ber in the Eversman radiation codes. Also, an equal sign has not been used here because in the actual 
Eversman codes a perturbation velocity potential mode amplitude, q>, is also included for each mode 
but does appear in Equation 34, Reference 3. Thus, including the perturbation velocity potential 
mode amplitude in the above equation: 

Vwnkpfrr't't) = * x ) (V.9) 
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The perturbation potential relates to the acoustic pressure through Equation 16 of Reference 3: 


PWnk Mr J X ' r ’ < t > ^ 



(V.10) 


where ^ is the modal acoustic pressure in the Eversman radiation codes. 

Substituting Equation (V.9) into Equation (V.10) results in Equation 122 in Reference 3. A form of 
this equation which includes the exponential is: 


Pw^J X ’ r ’M ~ ~ *>/" 


UfkxWn kfi 


0) 


(V.«) 


Non-dimensionalize this equation by r d and to get: 


PfC 00 

- - f — * 


1 - 


MfkxWnkn 


*Jr 


^Wnk^Wmfi^ 


- m<p k xWnkt p j 12) 


Note that the perturbation velocity potential mode amplitude, <p, the local static density, pp the duct 
radius, r d , and far -field acoustic speed, c TO are dimensional so that Pw„kn rad is dimensional. Compare 
this expression to Equation (20) which is the expression for the acoustic pressure in CUP3D. For a 
single mode Equation (20) is written as: 


P^r.M (V- 13 ) 


Comparing Equation (V.12) to Equation (V.13) it is observed that: 

• The axial origin of the Eversman Radiation codes may be different from that of CUP3D. 

• One equation appears to be the complex conjugate of the other, i.e. = Pwhkft 

Regarding the axial origin, in CUP3D the axial location of the origin is given explicitly as the interface 
plane so that at the interface plane, x - x p w = 0. In the Eversman radiation codes, the axial origin 
at the input plane is set to be at the interface plane, so thatx = = 0. As a result, the axial loca- 

tions are forced to be the same for CUP3D and the radiation codes. 

Now we check the complex conjugate observation by taking the complex conjugate of Equation (V.12): 



(V14) 


Note that the complex conjugate observation above is not true for cutoff modes. Instead, we must 
compare Equation (V.13) with Equation (V.14). These two equations cannot be considered the equal 
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to each other unless y Wnkfl - \kxw n knj • With some analysis (which will not be covered here), this 

expression has been found to be true for both the Eversman inlet and aft radiation codes. In the case 
of the inlet radiation code, a difference in the definition of positive* and Mf complicates this result 
slightly, but does not impact on the equivalence of the two equations. 

For consistency, make the expression in the exponential common between CUP3D and the radiation 
codes. Also, xp is a real number. Thus: 


PWnlv ~ ~ 1 ^ V> 


Mfh 


Wnkfi 


* 

\<PwJr)} (V.15) 


By setting the amplitudes of Equation (V.13) equal those of Equation (V.15) and solving for A P Wn)v ' 


.PfC< 

= l-r — 


Wnkfi p <x,r d 


Vr 


1 - 


MfkxWnkfi 


Vr 


4 


(V.16) 


Thus this expression for A p Wnk ^ is used to calculate the scattering coefficients for the radiation codes. 

V.2. 1. 1 Inlet Radiation Code Scattering Coefficients 

Referring to the scattering coefficients for the inlet radiation code in Section 1.3.2 of Appendix I, we 
note that only the Sff scattering coefficient is needed where D refers to the interface plane where the 
noise source and the inlet meet (e.g. in Figure 7 ,D = 1). 

Equation (22) shows how scattering coefficients may be calculated for the inlet radiation code. Thus, 
knowing that radiation codes only scatter on p, and assuming a unit amplitude pressure wave input 
(i.e.^ = 1.0 + i0.0): 


cDD 

21 nkfjL''/ikfi 


^2_ ^ 

A lnkp 


(V.17) 


Substituting Equation (V.16) into Equation (V.17): 




Vr 


4 


cDD 

°21 nku'/tkfi 



(V.18) 


Since this is a ratio, non-dimensional perturbation velocity potentials can be used in this equation. 
Equation (V.18) is the expression used to calculate scattering coefficients in the inlet radiation code. 
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V.2.1.2 Aft Radiation Code Scattering Coefficients 

Referring to the scattering coefficients for the aft radiation code in Section 1.3.5 of Appendix I, note 
that only the scattering coefficient is needed where U refers to the interface plane where the noise 
source and the nozzle meet (e.g. in Figure 7 ,U — 3). 

Equation (25) shows how scattering coefficients may be calculated for the aft radiation code. Thus, 
knowing that radiation codes only scatter on n, and assuming a unit amplitude pressure wave input 

a* - 10 + i00 > : 


cUU 

° 12 nkfi ' 'AkjJ. 


A 1 'nkf± 
^2 n kn 


(V.19) 


Substituting Equation (V.16) into Equation (V.19): 


cUU 

°12 nku'/ikfi 



(V.20) 


V.2.2 Radiation Pressure Directivities 

To correctly output far-field directivities from the Eversman codes, examine Equation 83 in Refer- 
ence 3. Putting this equation into the notation of this report: 

Vw+farM) " (V21) 

We can convert the coordinate system from cylindrical to spherical coordinates and change notation 
so that: 


The acoustic pressure for a given mode can be written in the same form: 

When comparing this equation to Equation (IV. 1) in Appendix IV we observe that: 

* 

**,<*.») - 

For a given wave type, Equations (IV.4) or (IV.5) give the directivity calculation. In the Eversman radi- 
ation codes where there is only one upstream going pressure wave and one downstream going pressure 
wave: 

Pnk/i(R’Q) = ^Wnl^PnkfSk^ (V.25) 


(V.22) 

(V.23) 
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So that in the radiation codes: 


® •* A w nkn m P ® ) (V.26) 

or: 

* * 

W* 9 > = J [*W*,*)] (V27) 

When the radiation codes are run, the input value is A*L , = 1.0 + t'0.0. Thus for this case: 

" n *t i rad. 

* 

( V28 ) 

Since is calculated by CUP3D after coupling, the radiation codes only need to output the direc- 

tivity shapes based on the above unit input. Therefore: 

* 

DIRECTIN(n,k,fi,IANGLE) = (V.29) 

where DIRECTIN may be found in Appendix I as an Acoustics Properties File parameter for inlet 
and aft radiation codes. 
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